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Abstract 


The  relative  significance  of  thermal  stresses  as  a las- 
er heating  damage  mechanism  is  assessed  by  comparison  with 
the  damage  mechSnisms  of  melting  and  thermal  degradation  of 
structural  strength  properties.  The  limiting  cases  of  one- 
dimensional axial  and  one-dimensicral  radial  heat  flux  in 
a thin  target,  plate  whose  plane  is  normal  to  the  axis  of  a 
stationary,  axially  symmetric  heat  source,  are  investigated. 
A on(  -dimensional  radial  heat  conduction  numerical  model  of 
the  linear  thermoelastic  stress  field  including  the  effects 
of  melting  and  structural  failure  is  developed.  Residual 
tensile  strength  and  damage  size  are  presented  as  functions 
of  the  laser  beam  and  target  plate  parameters. 


THERMAL  STRESSES  AS  A 


LASER  HEATING  DAMAGE  MECHANISM 


I.  Introduction 


The  advent  of  the  laser  with  its  inherent  potential 
capability  of  depositing  very  high  intensity  energy  on  a 
target  has  generated  interest  in  possible  applications  of 
laser  energy  to  cause  structural  deunage  to  a target.  This 
thesis  presents  the  results  of  a study  to  assess  the  signif- 
icance of  1' ser  induced  thermal  stresses  as  a damage  mech- 
anism in  a metal  target  by  compai;ison  with  the  damage  mech- 
anisms  of  melting  and  thermal  degradation  of  structural 
properties. 


Statement  of  the  Problem 

The  objective  of  .this  study  is  to  determine  whether 
Uiezinal  stresses  are  an  efficient  mechanism  for  causing 
structural  damage  to  metal  targets  relative  to  the  other 
laser  heating  induced  damac^e  sources  of  melting  and  thermal 
degradation  of  structural  properties.  The  study  is  anelyt- 
ical  and  based  primarily  orT  the  appi%c.ation  of  theoretical 
heat  conduction  and  linear  elasticity.  A major  goal  is  to 

" I . -,.»i  * j,  ■ 

develop  analysis  {hi’oeedof^s  tidtich  ere  ^relatively,  simple  to 
apply  (i.e. , Nlfioh  do ^t  requirt  t^a  usa  of  large,  expen- 
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sive  to  apply,  coa^uter  prograns)  so  that  the  procedures  may 
be  used  to  obtain  rapid  estimates  of  laser  heating  induced 
thermal  stress  fields  under  varying  conditions. 

Significance  of  the  Problem 

Most  of  the  studies  of  laser  heating  damage  conducted 
thus  far  have  been  experimental  and  concerned  with  target 
melting  and/or  material  property  degradation  (Ref  12-18,  25, 
30,  31).  No  known  unclassified  studies  have  been  made  to 
determine  quantitatively  the  relative  Importance  of  thermal 
stresses  as  a source  of  d2unaging  target  materials.  Thermal 
stresses  are  of  interest  because  of  the  possibility  that 
significant  structural  damage  can  be  induced  for  lower  laser 
beam  power  or  peak  intensities  than  for  the  other  two  damage 
sources  under  consideration.  A need  exists  for  an  analyt- 
ical procedure  to  rapidly  assess  the  relative  severity  of 
damage  to  be  Induced  by  thermal  stresses,  melting,  and 
thermal  degradation  in  metal  targets  as  functions  of  laser 
and  target  parameters.  Experiments  to  obtain  similar  data, 

* i 

particularly  concerning  thermal  stresses,  can  be  expected 
to  provide  only  minimal  data  in  the  near  future  because  of 
the  cost  and  time  required  to  conduct  sufficient  tests  over 
. V the  wide  range  of  parameters  of  interest,,  and  due  to  the 
di!!ff iculty . in  adequately  instrumeintiug  numerous  target  spec- 
imens to  obtain  accurate  tbermal  stress  data.  An  analytical 
« mode^  thus  could  serye  to  provide  such  information  in  lieu 
of  expeximents  and  as  a guide  t6  planning  efficient  experl- 
inents«in  the  futiXte. 
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Physics  of  Heating  Metals 
by  Laser  Radiation 

The  following  description  of  the  physical  process  by 
which  laser  light  radiation  is  converted  to  heat  energy  in 
a Betal  target  is  adapted  fron  Ref  1:2-7.  When  a laser 
beam  is  directed  at  a metal  target,  only  a fraction  of  the 
initial  beam  energy  is  absorbed  by  the  metal.  A portion  of 
the  initial  energy  which  is  in  the  form  of  photons  may  be 
attenuated  by  the  intervening  medium,  usually  air.  A large 
fraction  of  the  photons  incident  on  the  metal  target  may  be 
reflected.  For  the  problem  of  interest,  the  assumptions  are 
made  that  the  photons  absorbed  by  the  metal  are  converted  to 
heat  essentially  instantaneously  and  within  a very  thin  layer 
(relative  to  the  target  thickness)  at  the  surface  of  the 
metal.  These  assumptions  have  been  shown  to  be  valid  for 
aircraft  structural  metals  (aluminum,  steel,  titanium,  mag- 
nesium) irradiated  by  continuous  wave  lasers  (Ref  1:3-5). 

Following  the  conversion  of  the  light  energy  to  thermal 
energy,  heat  is  transferred  from  the  metal  surface  by  four 
mechemisms.  The  only  one  of  these  mechanisms  to  be  consider- 
ed«in  the  present  study  is  that  of  heat  conduction  within 
the  metal.  The  three  mechamiams  being  neglected  here  are 
radiation  from  the  metal  surface,  convection  in  the  adjacent 
medium,  and  removal  by  gravity  or  airflow  of  any  melted 
(liquid)  portion  of  the  metal.  Only  the  latter  mechanism 
(melt  removal)  affects  a comparison  of  the  relative  signif- 
icance of  the  damage  modes  of  interest.  The  efficiency  of 
the  melting  mechanism  obviously  is  dependent  on  the  behavior 
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of  the  melt  at  the  metal  surface.  For  this  study,  the 
assunq;>tion  is  made  that  any  melt  formed  is  removed  instan- 
taneously fay  an  unspecified  mechanism  so  that  the  most  ef- 
ficient mode  of  the  melting  mechanism  is  considered.  Studies 
of  melt  removal  are  contained  in  References  1,  22,  23,  25,  31. 

Scope  of  Study 

Limiting  cases  of  laser  .heating  induced  damage  were 
studied  by  considering  the  limiting  heat  conduction  cases 
of  one-dimensional  axial  (through-the- thickness)  and  one- 
dimensional radial  conduction.  The  study  is  theoretical, 
based  on  the  quasi-static,  uncoupled  thermoelasticity  theory. 
The  problem  is  that  of  a finite,  stationary  heat  source  nor- 
mally incident  at  the  center  of  a large,  thin  plate.  The 
plate  material  is  assumed  to  be  a homogeneous,  isotropic, 
linearly  elastic  material.  The  plate  material  properties 
are  assumed  to  be  temperature  independent  with  the  excep- 
tion of  the  strength  properties.  Rest  losses  to  all  sources 
are  neglected. 

For  the  axial  heat  conduction  problem,  analytical  solu- 
tions are  applied  to  obtain  temperature  diitributlons  for 
input  to  the  thermal  stress  problem.  Numerical  integration 
of  the  thermoelasticity  equations  is  used  to  obtain  the  ther- 
mal itress  distributions.  For  the  radial  heat  conduction 
case,  a numerical  model  of  the  heat  flux  based  on  the 
Fourier  heat  conduction  equation  la  applied. 

Milting  and  thermal  degradation  of  the  plate  material 
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structural  prqperties  are  inodeled  for  comparison  with  ther> 
nal  stresses  as  damage  MChanisms.  Damage  due  to  the  com- 
bined effects  of  these  mechanisms  is  also  considered. 
Measures  of  plate  damage  used  are  residual  tensile  strength 
and  materiel  failure; 
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II.  Outline  of  Study 


The  primary  objective  of  this  effort  was  to  detormine 
the  relative  efficiency  of  ihertvil  stresses  as  a damage  mech- 
anism due  to  heating  a target  material  with  a stationary  con- 
tinuous wave  (CH)  laser  beam,  nie  problem  basically  consists 
of  applying  the  concepts  of  both  heat  transfer  and  thermo- 
elasticity. In  the  basic  theory  for  describing » mathemati- 
cally, the  behavior  of  elastic,  isotropic  solid  media  under 
the  combined  action  of  heating  and  external  loading,  the 
problem  i-^  one  of  a coupled  thermal/mechanical  boundary 
value  problem  (Ref  5:3).  For  most  applications,  useful 
solutions  to  the  problem  require  that  certain  simplifica- 
tions be  made  to  the  coupled  theory.  The  usu^l  simplifi- 
cations result  in  eliminating  the  thermal/mechanical  cou- 
pling and  inertia  terms.  Elimination  of  the  coupling  be- 
tween mechanical  and  thermal  effects  enables  the  thermo- 
dynamic and  mechanical  parts  of  the  problem  to  be  analyzed 
separately.  If  the  inertia  effects  are  negligible,  the 
mechanical  part  of  the  problem  reduces  to  one  of  static 

i 

thermoelasticity. 

The  subject  study  began  with  an  investigation  to  deter- 
s:‘ne  a suitablp  form  of  the  thirmal/mechanioal  theory  for 
the  present  amplication,  as  discussed  in  Section  III.  Solu- 

A 

tions  to  the  thermal  and  mechanical  e<|uations  were  then 
investigated  to  describe  the  heit  transfer  and  thermal 
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8tr«8S  field8  for  the  8pecific  problem  at  hand,  lasi;r  beam 
heating  of  thin  platen.  The  modele  thu8  developod  are  dis* 
cussed  in  Sections  IV  and  V.  Sctction  VI  describes  the  models 
used  to  cdiaracterize  the  melting  and  thermal  degradation 
mechanisms  for  ccmparisons  with  the  thermal  stress  damage 
source.  Section  VII  contains  the  analysis  of  the  damage 
mechanisms,  input  parameter  sensitivity  and  relative  effi- 
ciencies based  on  the  one-dimonslonal  axial  and  radial  heat 
conduction  models.  Study  results  and  conclusiors  are  pre- 
sented in  Section  VIII. 

i^pendix  A is  a listing  of  the  one-dimensional  radial 
flux  niimerical  model  developed  as  part  of  the  i.tudy.  A 
sample  output  listing  is  given  following  the  program  list- 
ing. Appendix  B describes  the  algorithms  used  in  the  radial 
flux  numerical  code.  Appendix  C presents  the  :nalysis  con- 
ducted to  validate  the  radial  fliix  code. 

The  damage  effects  of  laser  heating  were  d ' .ermined  for 
the  ranges  of  laser  beam  and  target  parameters  lisced  below. 

1.  Absorbed  laser  beafn  powers  up  to  al:.ji't  30,000 
jottlGs/sec 

2.  Laser  beam  diameters  between  5 and  15  cm. 

3.  Target  plate  thicknesses  of  0.1  to  0.65  cm. 

4.  Large  plate  diameters  relative  to  beam  diameter 
A single  target  material  was  considered,  2024-T3  aluminum, 
with  a constant  absorptivity  of  0.5.  The  target  plate  was 
assumed  to  be  heated  in  Uie  absence  of' external  mechanical 

i 

loading. 


III.  Theory  Development 
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This  b&ckground  discussion  of  the  fundamental  mathcdnati- 
cal  theory  for  the  behavior  of  a homogeneous,  isotropic, 
elastic  solid  under  the  combined  action  of  heat  and  external 
loading  is  adapted  frcm  Ref  5.  This  behavior  is  uniquely 
describe!  by  the  following  four  equations  together  with 
proper  initial  and  boundary  conditions: 


w,ii 

« pCgTt  (3X'f2p)aT^£]^)( 

(1) 

®ij'j 

- p2l 

(2) 

(3) 

(4) 

where 

k » thermal  conductivity 

T » temperature  (absolute) 

T-  » reference  temperature  (absolute)  at  which 
material  is  stress  free 

p • material  density 

Cg  * specific  heat  at  constant  deformation 

vE 

1 » bame  constant  - 
M • L«ne  constant  6 
V M Poisson's, ratio 
fi  ■ Young's  modulus 
O G « shear  modulus 
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c » linear  strain 
a » stress 
u « displacement 

a » coefficient  of  thermal  expansion 


* Kronecker  delta 


!-  1 if  i*j 
* 0 if  if^j 


Indicial  notation  (lEKef  7:24-27)  is  used  here  for  brevity, 
however  familiarity  with  this  notation  is  not  necessary  for 
understanding  tlie  nuiterial  which  follows.  A brief  outline 
of  the  notation  is  given  here  to  clarify  the  subsequent 
discussion. 

Subscript  indices  are  used  to  refer  to  the  three  rec- 
tangular cartesian  coordinate  axes.  The  range  of  each  index 
is  therefore  three.  For  example,  the  three  coordinate  axes 
x^,  X2,  X3  can  be  expressed  as  where  i takes  on  the  values 
1,2,3.  An  indicial  equation  with  a single  index  on  each 
term  then  represents  three  equations.  A summation  convention 
is  employed  where  repeated  indices  in  a term  imply  summation 
over  the  range  of  the  repeated  index.  A comma  between  in- 
dices denotes  partial  differentiation  with  respect  to  one 
of  the  TOordinates,  i.e. 

*i,j 

Retuxming  to  the  subject  equations,  the  first  equation 
Is  the  energy  equation  for  the  linear  thermoelastic  theory. 
Equation  (2)  represents  the  equations  of  jwtion  which  reduce 
to  the  static  equilibrium  equations  of  elasticity  theory 
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«dien  the  inertia  tern,  pUj^,  is  zero.  Equation  (3)  is  the 
strain-displacement  relationship  and  equation  (4)  is  the 
stress-strain  or  constitutive  equation  for  an  isotropic 
solid. 

The  boundary  value  problem  governed  by  these  equations 
is  very  difficult  to  solve  in  general  (Ref  5:41).  For  most 
engineering  applications,  two  simplifying  assumptions  are 
made  which  reduce  the  problem  to  one  in  idiich  the  heat  con- 
duction and  thenooelasticity  equations  are  nor  coupled. 

This  allows  the  two  separate  problems  to  be  solved  indepen- 
dently. The  t«ro  sinq>lifications  are  that  the  thermal/me- 
chanical coupling  defined  by  equation  (1)  and  the  inertia 
term  (pu^)  in  equation  (2)  may  be  neglected.  Using  the 
terminology  in  Ref  5,  the  basic  theory  represented  by  equa- 
tions (1)  through  (4)  is  called  the  coupled  theory;  neglec- 
ting the  strain  rate  term  (£)^}  in  equation  (1)  results  in 
the  uncoupled  theory;  and  neglecting  both  the  strain  rate 
and  inertia  terms  produces  the  uncoupled  quasi-static  theory. 
•Whether  either  or  both  of  these  simplifications  are  appro- 
priate to  the  present  study  is  the  subject  of  the  remaining 

» 

discussion  in  this  section. 

A review  of  the  mathematical  problem  quickly  led  to 
the  conclusion  that  only  the  uncoupled  quasi-static  theory 
presented  reasonable  expectation  for  a successful  applica- 
ticm  to  the  subject  problem  within  the  temporal  constraints 
existing.  Hence,  the  studies  described  herein  are  based  on 

4. 

that  theory.  However,  the  implications  on  the  results  of 
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the  inertia  effects  nay  be  significant.  It  is 
shown  in  (Ref  5:43)  that  for  most  metals  (and  particularly 
for  aluminum  which  is  of  prime  interest  for  this  study)  the 
mechanical  coupling  tern  (strain  rate  term,  fj^)  in  equation 
(1)  is  negligible  if 

— « 60  (6) 
aT 

An  intuitive  argument  is  used  to  indicate  that  the  strain 
rate  should  bo  of  the  same  order  of  magnitude  as  the  time 
rate  of  change  of  the  temperature  if  there  are  no  sharp  vari 
ations  in  the  temperature  time  histories.  Since  the  strain 
rates  are  related  directly  to  displacement  time  histories, 
the  question  of  whether  mechanical  coupling  is  negligible 
is  related  to  the  magnitude  of  the  inertia  effects. 

It  is  further  concluded  (Ref  5:50)  that  the  mechanical 
coupling  and  Inertia  terms  are  negligible  if  the  rate  of 
heat  application  is  not  too  great.  Quantitative  ■* . Its  on 
the  application  of  the  uncoupled,  quasi-static  t)  - :y  must 
be  determined  for  each  problem  of  interest . Based  on  a 
review  of  example  analyses  (Ref  5:339,  4o€)  it  was  decided 

I 

that  such  an  analysis  wss  beyond  the  scope  of  the  present 

♦ 

study.  However,  it  I'A  probable  that  ior  at  least  the  higher 
heating  rates  and  short#;/  times  of  interfe.-«.r  i:>  tm  present 
study,  the  use  of  the  liucoupled  tn<v>v.v  is 

questionable  and  shoaid  her  subject  of  a sui/^<i:.;aent 
study  of  laser .-b-caa  thennal  stresses. 


■a 


IV.  Thermal  Models 


The  basic  situation  to  be  modeled  is  that  of  a station- 
ary laser  beam  normally  incident  on  a flat  plate.  The 
incident  beam  is  considered  to  be  axially  symmetric  so  that 
only  heat  flux  in  the  axial  and  radial  directions  occurs. 

That  is,  the  heat  'lux  is  assumed  to  be  two-dimensional  in 
cylindrical  coordinates  with  no  heat  transfer  in  the  circum- 
ferential direction.  The  effects  of  the  medium  between  the 
laror  beam  source  and  the  plate  are  not  included  as  these 
effects  do  not  contribute  to  the  physical  processes  of 
interest  for  this  study.  The  concern  of  this  study  beoins 
with  the  incidence  of  the  beam  energy  on  the  plate  <:arface. 

Attempts  were  made  to  obtain  a closed  form  solution  to 
the  problem  without  success.  Various  approximate  solutions 
are  given  in  Ref  4,  some  of  which  are  discussed  and  applied 
in  this  study.  Other  approximate  so?.utions  and  methods 
(Ref  26-29)  were  reviewed  but  all  involve  computational 
difficulties  or  complexities  which  are  contrary  to  the 
objective  cf  obtaining  solutions  that  are  relatively  simple 
to  apply. 

To  present  readily  tractable  problems  for  the  subsequent 
stress  analyses,  the  heat  transfer  problem  was  restricted  to 
modeling  the  two  limiting  cases  of  axial  and  radial  flux 
separately.  A one-dimensional  axial  flux  model  was  developed 
for  application  to  cases  of  relatively  short  heating  times 
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and/or  thick  plates.  The  other  limiting  case  associated 
^ with  long  hearing  times  and  thin  plates  is  thermally  modeled 

by  a one-dimensional  numerical  radial  heat  flux  model.  As 
discussed  in  Section  VII,  quantitative  definitions  of  these 
bounding  conditions  were  to  be  obtained  from  NASTRAN  (NASA 
Structural  Analysis) , (Ref  2)  a large  finite  element  program 
developed  by  the  National  Aeronautics  and  Space  Administra- 
tion (NASA) . This  goal  *:as  not  achieved  ccxnpletely  because 
o£  the  fairly  excessive  amcu!'t  of  time  required  to  run  all 
of  the  cases  necessary  to  determine  the  combinations  of 
heating  rates,  heating  times  and  plate  thicknesses  which 
could  be  modeled  accurately  by  either  an  axial  o:'  radial 

► 

flux  model. 

The  effects  of  heat  losses  from  the  heated  plate,  re- 

( ) 

' gardless  of  the  mechanism,  also  have  been  neglected  in  the 

present  study  in  that  they  are  not  expected  to  have  suffi- 
cient impact  on  the  results  to  justify  the  additional  com- 
plexity of  including  realistic  radiation  and  convection  loss 
conditions. 


Incident  Flux'  Models 

In  general,  the  spatial  distribution  of  the  incident 
beam  is  represented  as  a Gaussian  clstribution  which  is  gen- 
erally accepted  i\B  being  representative  in  analytical  studies 
of  laser  heating  effects  (Ref  3tl3,  75).  A uniformly  dis- 
tributed  beam  was  considered  in  the  axial  flux  cases  to 
simplify  the  st^.esb  analysis.  For  the  Gaussian  beam  the 


1 
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hMt  £lttx  density  absorbed  at  the  plate  surface  is  given  by 

-r*/2o* 

l(r)  * Ip^e  (7) 

idiere  Xp^  is  the  peak  absorbed  flux  density  at  the  beam  axis, 
r'  is  the  radial  distance  from  the  beam  axis  and  a is  the 
standard  deviation.  Xt-  is  convenient  to  define  the  beam 
radius  as 


a « 2a  (8) 

Therefore,  the  flux  density  becomes 


-2rVa* 

I(r)  « Ipj^e  (9) 

where  for  a'G&ussian  distribution  86.51  of  the  beam  energy 
is  contained  within  the  beam  diameter (2a) . The  total  absor- 
bed power  under  a Gaussian  beam  is  given  by 

w -r*/2o* 

« 2xlpj,  / re  dr  (10) 


^a  * 2irlp*®* 


(11) 


(12) 


The  portion  of  the  incident  beam  power  (P|^)  which  is  absor- 
bed at  the  plate  surface  is  given  ,by  the  absorptivity 


o 


0 


0 


Bence# 


pa 


20, Pi 


»a' 


(i4) 


and 


Kr) 


24, Pi  -2rVa* 


aa 


. * ^ 


(15) 


Heat  Conduction  Models 


The  purpose  of  the  heat  conduction  models  is  to  deter- 
mine temperature  distributions  to  which  thermal  stress,  mel- 
ting, and  thermal  degradation  models  may  be  applied.  Analyt- 
•leal  models  for  various  classes  of  axial  heat  Tlux  conditions 
were  investigated  and  a numerical  finite  thermal  element 
model  was  applied  for  the  radial  flux  case. 

Axial  Heat  Flux  Models.  Analytical  temperature  dis- 
tribution models  are  presented  as  solutions  to  the  differ- 
ential equation  of  heat  conduction  in  an  isotropic  solid 
(Ref  4:10) : 


pc 


8T 


ft 


9 

fx 


) * ly  ) * I? 


(16) 


where 


p **  material  density 
c ■ specific  heat 
T « t«i^rature 
t » tilM 

k •*  thermal  conductivity  > 

q ■ rate  of  heat  gain  per  unit  volume 


ir 


For  honogeneous  oolidt  such  that  tharaal  conductivity  does 
not  vary  with  poaiticMi,  tha  heat  conduction  aquation  re- 
duces to 


icv*T  SL  * 35 
pc 


idiere 


ic  ■ k/pc  * thermal  diffuaivity 
V " Z«aplacian  operator 


V*T  « + iil  + 

3x*  3y*  3*^ 


Isotr6pic«  homogeneous  solids  wore  the.;  only  type  considered 

in  this  study.  . An  additional  restri/;tion  applied  throughout 

^ * 

the  heat  transfer  'analyses  is  that  the  material  heat  trans- 
fer  proj^rties  are  not  fjonctions  of  temperature.  The  jus- 
tification for  this  restriction  is  that  the  complexities 

•»  ij 

introduced  into the  analyses  by  including  variable  thermal 
properties  were  not  believed  to.  be  necessary  for  this 
initial  st\idy« 

Closed  form  analytical  solutions  for  several  one-dimen- 
sional axial  heat  flux  cases  \mre  reviewed  for  possible  ap- 
plication. For  the  case  of  a semi-infinite  solid  with  uni- 
form flux  into  one  surface  and  insulated  on  the  opposite 
surface,  as  shown  in  the  following  sketch,  the  taiiperature 
distrlbuticS)  is  given  (Sef  4 till) 


C' 


I 


T(z,t)  * Tq 


^ ♦•fc  ♦.{32*  - l») 

+ ^ + SH 


k2* 


-ictn*ir*/t* 

e 


(19) 


0 


where 

♦•  *■  absorbed  heat  flux  density 

Tq  initial  uniform  temperature  throughout  slab 

For  thermally  thin  plates  or  long  enough  heating  times  such 
that  Kt/i*  > 1,  the  series  terra  becomes  negligibly  small 
and  the  thermal  model  reduces  to 

T(z,t)  - Tq  - 32*  - 1*)  (20) 

^Implications  of  those  models  to  the  general  problem  of  laser 
heating  of  solid  plates  is  discussed  in  Section  VII. 

Radial  Beat  Flux  ttodel.  For  relatively  low  heating 
rates  and  thin  target  plates  the  heat  flux  might  be  expected 
to  be^xne  eseentially  radial  after  acme  initial  period  in 

which  the  flux  is  primarily  exial  and  ths  plate  volume  under 

' - ^ 

» * 


I 


1 

the  bean  Leconea  heated  to  an  essentially  uniform  tcmqpera*  { 

] 

ture  through-the- thickness « One*dimensional  radial  flux  | 

I 

develops  a plane  stress  state  tihich  is  particularly  amenable  I 

to  analysis,  as  discussed  in  Section  V.  For  this  case,  a I 

numerical  finite  thermal  element  siodel  was  developed.  The  | 

1 

numerical  SK>del  also  allows  the  Gaussian  beam  profile  to  be 

! 

considered.  Figure  1 is  a schematic  illustration  of  the 
ntimerical  model  which  is  obtained  by  dividing  the  plate  into 

j 

concentric  cylindrical  elements.  | 

1 

The  numerical  model  is  used  to  solve  the  Fourier  heat  I 

conduction  equation  for  each  el«nent:  | 

‘ « 

I 

i 

-Jc  - -m  <w>  j 

• I 

where  I 


«>  rate  of  heat  flux  per  unit  area 

Heat  is  transferred  between  any  two  points  in  a solid  body 
only  by  conduction  and  equation  (21)  is  the  law  of  heat  con- 
duction for  isotropic  bodies  (Ref  5tl37).  A finite  differ- 
anoa  version  of  equation  (21)  is  applied  to  the  heat  flow 
between  adjacent  concentric  elements  as  indicated  in  the 
following  sketch} 


1) 


X Z-fl 
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ylhmre  Q • « total  flux  into  eleiBent  (I  -f  1) 

Banco,  applying  equation  (21)  at  the  interface  between  two 
elemente 

Q|?  ♦ 1)  . ^k(l  + 1)^4  k(l)]  [T(D  - T(I  + D] 

Ar(I  »1^  -t-  Ar(I) 

where 

A(Z  4-  1)  » 2tR£ 

R • inner  radius  of  element  (I  -f  1) 
i ■ plate  thickness 

Ihe  conductivity  and  incremental  radius  terms  in  brackets 
represent  averages  of  the  respective  quantities  taken  be- 
tween the  two  elements. 

For  the  present  study,  the  concentric  cylinders  all 
have  the  same  width  (Ar  » constant)  &nd  tlie  thexnnal  con- 
ductivity is  assumed  to  be  independent  of  position  and 
temperature*  Thus  equation  (22)  reduces  to 

Q(I  4 1)  - k(T(I)  - T(I  + l)  J2irRVAr  (23) 

or,  in  FORTRAN  code 

QXN(X  4 1)  « COND*(TBMF(X) 

-TBW^d  4 1))«2*FX^^£{./DELR  (24) 

The  tes^rature  of  each  element  ie  determined  by  itera- 


(22) 


tion  over  «guilX  tlae  increnento  of  an  energy  balance  equa- 
tion for  each  eleawnt,  at  indicated  in  the  following  sketch. 


OMlSd) 


QIN(l) 


QOOT(I) 


where 


OhBS(l)  ■ heat  energy  absorbed  into  element  (I) 
per  unit  tine 

QOOT(I)  « OIM(I  +1) 


Hence r the  net  increase  in  heat  energy  per  unit  time  is 

' OhBS(I)  -f  OXlUX)  - QIN(I  ^ 1) 

« 

This  increase  in  heat  .energy  results  in  a corresponding 
increase  in  the  internal  heat  ei\ergy  of  the  element.  As- 
susaing  conatant  density  and  specific  hetft,  i^e  incremental 
ii^rease  in  internal  heat  energy  per  unit  volume  is  given 
Henoev  the  energy  balance  gives  , , 

OhBS(z)  * QiNU)  -"omd  + 1) 

• ' ♦ 

f * 

» Rli0*C#*V01i*D^«KMPd)A>l^T  ^ ' (^5] 


Where 


hHO  » p » density 

■/*  ^ 

CP  • o « speoifio  beet 


, .1- 


i -s 


OiLTBfPCI)  » AT  ■ change  of  tmperature  In  tine 

incr«Mnt  At 

DELT  ■ At  « tine  increiaent 

VOL  « volmie  of  eleaent  (I)  « PI*EL* 

(2*1  - l)*0£LR*f2 

DSLTEMP  is  the  incrcnnental  tosperature  change  in  element  (I). 
This  may  be  expressed  as 

OBLTEMP(I)  - TEMP {1)2  - TEMP (I),  (26) 

* y 

nhere 

TEMP  (I)  ^ - oirrent  temperature  of  eluent  it) 

* 

TEMP (I) 1 « previous  t^operature  of  element  (I) 

« 

HencOf  equation  (25)  may  be  solved  for  the  current  tempera- 

* *“ 

turd: 

T^(I),  « TEMP(I)i  + DELT*(QIN(I)  + QABS(I) 

-QIN(I  + 1))/(CP#RH0*PI*BL*(2*I  - 1)* 

" DELR**2)  » V (27) 

I 

Equations  (24)  and  (27)  are  solved  for  each  thenhal  element 
at  each  time  increment "to  give  thp  radial  flux  temperature 
distribution.'  This  model  is  Subroutine  TBNPRAO  in  the 
TSTtmsS  progruR  ^sted  in  Appdnd|.x  A. 

I 

M maintain  heat  flow  iM:ability  in  such  k nun»ricalr 
incremental  time  model,  the  time  increment  must  be  kept  be- 
llow a certain  limit.  It  the  itime  In^nmaent  becomes  l#rge 

enough,  sufficient  heat  would  fl^  betweed  adjacent  elements 

•V- 

bo  cause  the  toi^rature  of  the  ini^ally  <K>oler  elmsent  to 


o 


1 

tmcome  larger  than  tha  initially  hotter  elteient  causing  the 
heat  flow  to  becone  reversed.  To  Maintain  heat  flow  stabil- 
ity the  cos^tational  tiae  incresant  for  the  one-diaiensional  ‘ 

radial  flux  Model  Must  be  restricted  to  (Ref  6:12) 

At  < ||1  (28) 

As  a chec)c  and  to  be  slightly  conservative,  the  TSTRSSS 
code  limits  the  tine  increment  to  ninety  percent  of  the 
value  defined  in  egn  (28). 

0 ! 


■» 


> 


0 


23  ■ 
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V.  Thermal  Stress  Models 


In  the  uncoupled  quasi-static  theory  the  thermal  stress 
prdi>lem  is  solved  in  sequence  after  the  heat  conduction 
problem  is  solved  for  the  tenperatvxe  distributions.  The 
thermal  stress  problem  is  uniquely  defined  by  the  static 
thermoelasticity  equations 

® 

c^j  * l/2(u^,j  + Uj,^)  (30) 

Oij  - + 2ycj^j  - (3X  + 2v)(x6^^{r  - T^)  (31) 

along  with  proper  boundary -conditions*.  For  the  case  when 
the  heat  flow  in  the  plate  is  primarily  axial  (i.e.  through- 
thc-thickness) , the  only  portion  of  the  plate  being  heated 
significantly  is  that  directly  under  the  beam.  For  a uni- 
formly distributed  beam  the  temperature  variation  will  be 
through~the**thicknes8  only.  Near  the  beam  axis  the  thermal 
stress  field  eah  be  approximated  by  that  for  a semi-infinite 
plate  of  thickness  2hf 


2h 


24 


Sinci^  4re  no  surface  tractions  on  the  edges  of  the 

pjletc  and  the  te??^ratur.«  la  a ftmction  of  z only,  and  due 
to  axial  3yrr;.!*etry,  it  is  expected  that 

®jj3£  * ®yy  * Hx)  (32} 

/ 

at  any  ti{*«.  Since  the  to^rj^ture  is  uniform  in  the  x - y 
plane,  no  shear  streases  are  generated: 

®xz  * ®yx  * ®zy  * ® (33) 

Also  since  no  constraints  on  expansion  in  the  z - direction 
exist, 

Czz  • 0 (34) 

All  of  the  bovndary  conditions  are  seen  to  be  in  terms  of 

r 

stresses  only.  For  thermoelasticity  problems  of  this  type, 

the  solution  is  facilitated  by  reformulating  the  thermo- 

elasticity  equations  in  temw  of  str^f^es  alone.  That  is, 

tha  displao«Bants  and  strains  ara  aliminated.  Iliia  is  dona 

* ^ 

thro&gh  the  use  of  the  strain  compatibility  equations  (Ref  7: 
124 } which  are  mathematical  constraints  that  insure  the  in- 
tegrabillty  of  aquations  (30)  to  c^tain  the  displacements. 

Xn  -Jis  stress  formnlatioii  of  the  problem,  the  equilib- 
ritm  equations  (29)  are  onphsnged^  but  stress-strain  aqua- 
tiont  (31)  and  ^ strain-cooipat^bility  aquations  are  re- 

f 


placed  by  the  so-called  stress  coapatibility  equations  (Ref 
5;d9): 

(1  + + as 

+ T,ij]  = 0 (35) 

For  th.e  present  problem,  the  equilibrium  equations  (29)  are 
identically  satisfiec  and  equations  (35)  are  satisfied  if 


Integration  of  (36)  yields 


■(S^) 


T + c,  + c.z 
1 2 


where  it  is  recalled  that  T » T(z,t)  - T^.  To  obtain  a nvon- 
trivial  solution,  Saint  Venant's  principle  (Ref  5:270)  may 

4 

be  applied.  That  is,  tha  requirement  for  zero  surface  trac- 
tions  at  the -edges  of  the  plate  is  replaced  by  the  statically 
equivalent  requirement  that  the  net  force  and  moment  pro- 
duced  by  some  a^x  distribution  oa  the  edge  he  zero.  These 
boundary  conditions  expressed  for  net  force  and  moment  per 
> unit  length  ef  the  plate  edge  bectxoe,  respectively 


is6 


(39) 


/ 


xds  » 0 


Substituting  (37)  into  (38)  and  (39)  gives 


Hence,  from  equation  (37) 


a 


XX 


/ Tdz 
h 


+ /**Tzd2 

^ -h 


(40) 


(41) 


(42) 


The  appl.vcfation  of  equation  (42)  to  the  present  study  along 
vith  suitable  thermal  models  is  discussed  in  Section  VII. 
This  result  is  given  in  (Ref  5:278). 

The  other  thermal  stress  models  of  interest  are  as- 
^ socinted  ./ith  conditions  when  the  temperature  distributions 
in  the  heated  plate  are  primarily  radial  only,  i.e.  constant 
throUgh-'the-thickness.  For  a thin  cylindrical  plate  with 
with  diameter  >>  thiclcness  and  no  surface  tractions,  the 
problem  is  one  of  plane  stress  such  that  the  only  non-zero 
stresses  ara 

dyj.  **  f I (r)  ^ (43) 

000  • f»(r) 

^7 


(44) 


using  cylindrical  coordinates. 

The  only  eguilibriun  equaticm  (Ref  5t248)  vdiich  £a  not 
identically  sero  is 


*®rr  . ®rr  "■  *^06 
“5F‘  * r • ® 

(45) 

From  the  stress-strain  eq^tions  (Ref  5:245) 

• 

Eejfj.  * a^j,  - VO00  + EoT 

• 

(46) 

“ ®00  “ ^rx  * 

(47) 

and  from  the  strain-displacement  equations 

> 

« la 

*^rr 

r 

(48) 

€09  - ^ (49) 

V 

Elindnating  Cga  between  (46)  and  (47)  and  substituting  (48) 

, * 

and  (49)  gives 

^ i v»)  - jlHlj  (50) 

Similarly,  eliminating  in  (46)  and  (47)  gives 

[•S'f  -oT(l,+  v)j  “ (51) 

Substituting  (50)  and  (51f  into  (45)  results  in 


fi  . (X  +v), 


Integrating  and  solving  for  the  radial  displacement 


u . 12022“  /'irdr  + c,  5 + ^ 
r a * I r 


Substituting  back  into  (50)  gives 


— f-  ^ Trdr  + (1  +v)~ 

l-\r  L r*  « 2 


- - '■>3 


By  substituting  the  boundary  conditions  of  zero  "<trac~ 
tions  (a^r  * 0)  on  the  inner  (r  « a)  and  outer  (r  «=  b)  radii 
of  the  cylinder r the  constants,  Ci,  and  cz  are  determined. 
Hence,  (Ref  5s290), 


“I  ri^-ail  Trdr  - /’ 

l.(b*  - a*)  a a 


Trdr 


oE  fir*  -f  a^) 
?■  L(b*  - a*) 


*>  r n 

/ Trdr  / Trdr  - Tr* 
a a J(56) 


For  the  particular  case  of  a jiolid  plate,  a •>  0 and 


[1  b 1 > 

ij.  / Trdr  - p-  ^ Trdr 


^Trdr  + i- 

► r* 


^*^Trdr  - tJ 


Another  particular  case  of  interest  is  a solid  cylin* 


Oi 


drical  plate  of  ehich  only  a finite  J-nner  cylindrical  por> 
tion  if  heated  appreciably.  If  it  is  assumed  that  this 
tested  core  (r  * aj)  is  at  a uniform  tempex^ture  (T<«  Tt) 

with  T » 0 for  r > ai,  thesl  from'  (57)  and  (58)  respectively, 

« 

f 

for  r £ a|, 


®06| 


oETi 


- 1 


and  for  r >,  ai. 


2 \b*  rV 


cETiaf 

2 


The  stresses  in  the  heated  core  are  seen  to  be  compressive. 
The  peak  tensile  stress  occjirs  jtSst  Outside  the  hOatad  qore 

and  is 


r « a, 


Ite  oaacitftA  shear  stress  ^ 

axliil  stress  fielf^  occurs  on  planes  at  i 45  degrees  to  the 

k " ' 

principal  stress  axes  and  is 


®rri 

T Ift  --  - * 

^Ultx 


in  the  heated  core.  Outeide  <^e  heated  core  the  maximum 
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shear  stress  is  given  by 


(64) 


Applications  of  these  thermal  stress  models  to  the 
laser  heating  problem  are  disctissed  in  Section  VI 1. 

fc. 
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VI.  Welting  and  Thermal  Degradation  Models 


To  assets  the  relative  effects  of  thexnal  stresses  as 
daiiiage  mechanisms  froei  laser  heating  of  metal  targets » the 
melting  and  thens'il  degradation  damage  mechanlmis  were  se- 
lected for  purposes  of  cos^rlson.  Other  damage  modes  such 
as  explosive  boiling » explosive  vaporization,  and  thermo- 
elastic  shock  waves  (Ref  12)  occur  at  generally  higher 
heating  rates  than  the  mechanisms  of  interest  In  this  study. 
Since  the  concern  Is  with  determining  efficiencies  of  the 
damage  mechanisms,  only  those  believed  to  occur  at  the 
relatively  lower  bean  powers  were  considered. 

Welting.  The  analytical  model  used  to  predict  melting 
effects  follows  that  of  (Ref  8:18).  Melting  is  assunw'd  to 
occur  over  a range  of  tonperatures , Ti  to  Tu,  where  is 
the  temperature  below  which  the  material  behaves  as  a solid 
and  T^  is  the  temperature  above  which  the  material  behaves 
as  a liquid.  For  the  present  study,  the  melted  material  is 
assumed  to  be  removed  innediately  upon  reaching  so  that 
the  effects  of  melt  retention  are  not  considered. 

During  the  phase  change  frcxn  solid  to  liquid  an  ef- 
fective specific  heat,  Cp«  Is  assumed  such  that  the  heat  of 
fusion  over  the  melting  range  is 

®pe  * ®p.  "*  for  Tj  < T Ty  (65) 


«duire 


Cp  “ specific;  heat  of  solid  pl.ase  (T  < Ti) 

£ib  **  heat  of  fusicm 

Once  any  portion  of  the  heated  target  reaches  the  upper 
Bttltihg  point,  that  portion  of  the  target  is  assumed  to  be 
instantaneously  removed  and  is  no  longer  considered  in  sub~ 
sequent  ccxiqmtatione . Damage  due  to  melting  is  then  defined 
the  reduction  in  target  strength  due  to  the  reduced  load 
carrying  area. 

Thermal  Degradation.  The  definition  of  thermal  degrada> 
tion  as  ai^lied  in  this  study  is  the  reduction  in  structural 
strength  properties  with  increasing  temperature  of  thin 
plates  as  described  in  Ref  9.  Although,  as  discussed  in 
Section  VIZ,  the  primary  measure  of  damage  being  used  for 
this  study  is  the  residual  tensile  strength  in  a uniaxial 
loaded  panel,  the  effects  of  thermal  degradation  on  tensile, 
CK>n^res8lve  and  shear  strength  properties  are  considered. 

The  only  material  considered  is  2024-T3  aluminum.  Ref 
9 contains  data  for  the  effects  of  temperature  on  material 
strength  as  a fuwation  of  time  at  temperature.  The  parti- 
cular data  used  here  is  contained  in  Figures  3.2.3.1.1(a)., 
3.2.3.1.2(a).,  and  3.2.3.1.2(b).  of  Ref  9 for  ultimate  ten- 
sile strength,  con^ressive  yield  strength,  and  ultimate 
shear  strength,  respectively.  These  figures  are  reproduced 
here  as  Figures  2,  3,  end  4. 

Zt  was  decided  to  use  the  helf-hour  exposure  data  in 
itef  9 after  reviewing  available  data  on  the  more  rapid  heat- 
ing effects  associated  with  high  i^nrer  laser  heating.  Such 
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Pig*  3 Thermal  Degradation  of  Compressive  Yield  Strength 
for  2024-T3  AltMinum 
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Fig*  4 Thermal  Degradation  of  Ultimate  Shear  Strength 
for  2024-T3  Alunimai 
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date  are  contaii^  in  Ref  13  through  18.  Doe  to  dai:a  scat* 
ter,  variaticms  in  teat  nethods,  inconsistent  teat  pari^* 
■aters,  etc.,  these  data  do  not  pres«it  a clearly  defined 
relationship  bettfeen  stxength  and  exposure  time.  Additional- 
ly the  data  for  the  aaterial  of  interest,  2024-T3  aluainuis, 
are  concerned  only  with  tenaila  properties.  For  the  present 
study,  degradation  of  cosg>ressi<»i  and  shear  strengths  were 
also  of  interest.  The  theneal  data  of  Ref  9 provides  the 

9 

necessary  tensile,  coaqpressive  and  shear  degradation  rela* 
ticmships  for  consistent  conditions.  Although  the  data  for 
the  much  shorter  heating  tioies  presented  in  Ref  13  through 
18  indicate  that  strength  properties  for  the  short  heating 
tines  associated  with  laser  heating  can  vary  appreciably 
fron  the  half-hour  exposure  data,  it  w^s  decid^  to  use  the 
latter  data  for  the  present  study  based  on  the  presuniption 
that  these  data  are  suitable  first- approximations  and  they 
are  of  standardized  validity,  complete,  and  consistent 
which  at  least  partially  c<xnpensate8  for  possible  detailed 
discrepencies  with  the  shorter  heating  time  effects. 

Based  on  the  half-hour  exposure  curve  of  Fig.  2,  for 
cooq^utational  convenience  it  was  decided  to  approximate 
the  ultimate  tensile  strength  by  two  straight  lines  such 


that 


*^ttt  • Fto  * C < T < 200*C 


(86) 


" 0.75Fto  - 200(T-200) 


300*C  < T < 370*C  (67) 


36 


'where 


« 0 T > 370*C 


« oltlMite  tpusile  strength  at  teoperature 
T(n/aii*) 

Pi^  « tkltisMte  tens  ile  strength  at  roon  teaqperature 
- 44»820  n/csi*  lor  2024*T3 

¥ « tesperatntWf'  *0  above  room  tei^>erature 

Similarly,  the  half-hour  ejqpostire  compressive  yield  strength 
was  modeled  as 


'co  r 25T 


0 < T < 232»C 


- 0.75Fqo  - 140<T-232) 


232*C  < T < 370*C 


T > 370*C 


vrhere 


4 h 


Pey  ^ conqpresslve  ‘yield  strength  at  temperature 

■ room  twqperature  ccMipi^ssive  yield  strength 
•»  23,440  n/cai*  for  2024-T3 


•She  ultimate  shear  *htrength  was  mpdeled  analyticaily  as 


S' 

F»a  • 0 < T < 216^C 


Fgu  • d*75Fg^  - U2(i-StX6)  ' 2X6*C  < T < 370»C  (73) 
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* 0 T > 370*C  (74) 

6 


• ultiaatd  shttar  strength  at  temperature 
T(n/c«*) 

^so  **  ttltiaate  shear  strength  at  roon  teaqperature 

• 25,510  n/qs*  for  2024-Y3 


VIX.  La«.ir  B— ting  Effects 


f 


V 

One»Dlaen«lonal  Heat  Flux  Model  Limits 

Since  the  pi^sent  rtudy  was  concerned  with  the  Uniting 
cases  of  axial  and  radial  heat  flux  conditions,  an  initial 
question  to  be  answered  was  under  what  conditions  were  the 
assuoq»tions  of  <me»dinertsional  axial  or  radial  flux  valid. 

It  was  originally  planned  that  the  ranges  of  input  para> 
neters  for  which  the  heat  flow  could  be  considered  to  be 
essentially  one-dinensional  (either  axial  or  radial)  were 
to  be  detemined  by  coBg>aring  temperature  and  stress  dis- 
tributions fron  the  one-dinensional  models  ' those  obtained 
from  a large  two  or  three-Alssnsional  analysis  finite 
element  cenqiuter  program. 

Initially,,  the  MABC  (Ref  10)  three-dimunsional  trans- 
ient thermal  stres^  code  was  to  be  used  but  was  abandoned 
when  access  to  the  program  was  terminated  by  non-renewal  of 
a lease  contract  between  the  Air  Force  Plight  Dynamics 
Laboratory  and  Control  Data  Corporation.  Efforts  were  then 
made  to  apply  the  NASTRAN  (Ref  2)  program  using  the  Navy 
thermostructural  analysis  additions  and  modifications  (Ref 
11).  C<msiderable  time  was  expended  in  getting  this  program 
t6  function  properly  due  to  laoh  of  familiarity  with  NASTRAN 
•nd  doa  to  insufficiant  documentation  in  Ref  11*  As  a con- 
sequence of  these  anA  other  difficulties  with  the  program, 
it  was  daoidsd  to  aliminate  this  approaoh  to  dafinlng  tha 


-one^ijisnsioiMlity  liaits.  As  discasssd  later  in  thifi  sec- 
RASTRAM  was  used  to  vallate  the  niueerical  one- 
diaensional  ri^ial  flux  cos^tations- for  tenperatures  and 
themal  stresses. 

An  alternate  approach  was  used  to  define  the  limits  of 
the  one-dimensional  axial  and  radial  flux  thermal  stress 
models  based  on  the  results  of  the  studies  presented  in 
Ref  8 and  19.  There  is  reported  the  development  of  a 
single  parameter  which  is  shown  to  determine  the  limits  on 
the  one-dimensionality  of  heat  flux  for  the  stationary  beam 
melting  problwa.  The  parameter  is  a dimensionless  absorbed 
power  per  unit  thickness  defined  by 


^la 


i^PrCp^Tm-To)4.Lml 


(75) 


idiere 


I » plate  thickness 
^m  a melting  temperature  ° 

Ljq  a heat  of  fusion 

In  Ref  19  it  is  shown  tl^at  i^elt-throu^  times  based  on  the 
assunption  of  one^dimeniional  axial  conduction  closely  pre- 
dict those  ^roit  a two-oilmen sional  numerical  analysis  for 
values  of  <^reater  jbhan  about  70.  Values  of  less 
t^  about  5 sre  siiown  td  indicatm  that  the  flux  ia  primari< 
ly  radial.  In  (Ref  30)  an  oeffeetiva*  dimansionless  power 
per  unit  thic)cneas  ii  developed  for  correlation  with  melt- 
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throngh  tiae«  which  ii|clude8  an  i^roxinate  allowance  for 
the  t«ag>eratiire  dependence  of  Cp. 

Although  these  one*^ijKmsionality  liaits  are  based  on 
the  Melting  prcd>lesi»  it  is  reaso^ied  here  that  the  assun^tion 
that  these  liaits  would  also  apply  for  the  thermal  stress 
pr^l«B,  at  least  for  heating  times  up  to  melting,  is  con- 
sistent with  the  earlier  assus^ion  (Section  III)  that  the 
strain  rate  is*  of  the  same  order  of  magnitude  as  the  ti':ie 
rate  of  ch2mge  of  t^aperature.  That  Is,  it  seems  reasonable 
to  assume  that  the  one-dijeaensionality  limits  for  the  heat 
^<mduction  problem  (based  on  melting  times)  should  be  about 
the  same  as  for  the  thermal  stress  problem  if  the  strain 
rates  closely  follow  .the  tespereture  tiae  rates.  Conse- 
quently,  the  above  noted  limits  on  one-dimensio \al  axial 
emd  radial  flux  are  applied  in  the  present  study. 

For  the  material  considered  in  this  study,  3024-T3 
aluminum,  the  one-dimensionality  limits  based  on  can  be 

i 

reduced  to  limits  on  laser  beam  po(lrer  and  target  thickness. 
For  2024-T3  aiuminiUR,  r 

i * 

r *“ 

K ■ 0.514  cm? /sec  , 
p ••  2.78  |m/cm*  > 

Cp  • 1.05  joulp/gm-*C  . 

Tj,  - 

To  • 24*C 

« 375  5oule/gm 

^»a  " Ijffr  <”* 

41. 


f 


C) 


Hence » for  cme-diMeneional  exial  conduction  ^ 

^ > 8.55  X 10'  (axial  flux)  (77)  ] 

Jw  I 

i 

Similarly,  for  one>dimensional  radial  conduction,  | 

< 5 and 

^ < 6.11  X 10’  (radial  flux)  (78) 

r - ‘ I 

Axial  Plux  Analyes 

For  relatively  high  heating  rates,  the  initial  heating 
of  a plate  with  the  beam  incidence  normal  to  the  plate  mid- 
plane would  be  expected  to  be  primarily  axial  through-the- 
thicicness  of  the  plate.  Since  the  only  area  being  heated 
significantly  is  that  directly  under  the  beam,  the  thermal 
model  of  equation  (19)  gives  the  temperature  distribution 
through-the-thicKness.  To  obtain  a closed  form  analytical 
solution  to  the  thermal  stress  problem,  the  approximation 
for  thermally  thin  plates  (Kt/t’>l)  is  applied  so  that  the 
temperature  distribution  is  given  hj,*  equation  (20).  If  this  j 

thermal  model  is  applied  to  the  thermal  stress  model  of  i 

! 

equation  (42) , the  in-plane  normal  stresses  near  the  beam  j 

axis  are  given  by 


®xx  • ®vv  “ ■ — (h*-3s*) 

**  12(l-v))ch 


(79) 


O 


where  the  coordinate  systiwn  of  Section  V is  used.  This  sol- 
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tttioR  is  se«n  to  be  independent  of  tine  once  the  initial  re- 
striction that  Kt/A^  > 1 is  set.  Hie  naxianm  tensile  stress 
occurs  at  s » 0.  The  stress  distribution  through  the  plate 
is  shown  in  Fig.  i 

An  intenRediace  thermal  situation  can  be  hypothesized 
%diich  might  approximate  the  peak  stresses  during  transition 
from  a predauinantly  axial  flux  case  to  ohe  in  %diich  the 
flux  becomes  essentially  radial.  After  the  core  area  under 
the  beam  is  heated  for  sor  time  by  axial  flux,  this  heated 
core  will  become  heated  to  an  essentially  uniform  tempera- 
ture. At  this  time  and  prior  to  the  development  of  signi- 
ficant radial  flux,  the  thenaal  stress  problem  can  be  con- 
sidered as  one  pt  a heated  center  disk  surrounded  by  an  un- 
heated elastic  cylinder.  Equations  (59)  through  (64)  thus 
provide  the  thermal  stress  model  for  this  problem.  An  ini- 
tial assessment  of  the  stresses  generated  for  this  case  is 
obtained  from  equation  (62)  for  the  peak  tensile  stress 


(80) 


Ti  is  the  \miform  temperature  of  the  heated  core  as  develop- 
ed  from  equation  (20) , 


T(z,t)  - « iiU.i6Kt+3z*-A*l,  kt/A*  > 1 (81) 

® m 


: 


0 


Arbitrarily,  to  define  the  heating  time  when  the  flux  vec- 
tor  would  be  expected  to  be  primarily  radial  let  it  be 


^3 


I 


required  that  the  front-to-back  eurfaee  tei^rature  varia> 
tlon  be  restricted  to  10%  or  less  of  the  frmit  surface  (s  « 
i)  teaperature.  Froa  equation  (81)  above 


#t)  - To  - ^Oict-H*) 


The  bac)c  surface  (s«0)  te^^ratura  is 


T(0,t)  - Tq  • 


Subtracting  (83)  fron  (82)  gives  the  frcmt-to-back  ten^ra- 
ture  variation  f iTfi.B 


Requiring  that 


results  in 


i 0.l[TU,t)-ToJ 


< -iiL(3Kt+i*) 

2k  - 30kt 


which  reduces  to 
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This  r«8ttlt  it  teen  to  be  an  increate  in  the  themal  nodel 


rettriction  of  ict/i'  > 1. 

The  toltttiont  for  axial  flux  probleet  given  that  far 
ai^ly  to  l<mg  heating  tinea  and/or  thin  platet  (Kt/i*  > 1). 
POr  earlier  heating  tinea  (or  thicker  platet)  the  thermal 
node!  of  eqn  (19)  givet  the  temperature  diatribution  through- 
the-thickneaa  of  the  plate.  The  thermal  atresaea  can  then 
be  determined  by  numerical  integration  of  the  strata  model 

X 

given  aa  eqn  (42).  Pig.  6 tuanarites  the  resulta  of  such 
an  analytia.  It  waa  initially  fotmd  empirically  that  the 
in-plane  atreaa  it  given  uniquely  by  the  value  of  <t/t* 
at  any  depth  ix/t)  in  the  plate  aa  ahoim  in  Fig.  6.  This 
eaa  aubaequently  c^firmed  analytically  by  factoring  (4«t) 
out  of  the  thermal  model  equation  giving 


(88) 


» 4 1 

The  combined  factor  li-  has  the  dimensions  of  temperature 

k 

and  the  terms  inaide  the  bracketa  are  dimenaionleaa.  Exam- 
ination of  the  thermal  stress  equation  (42) 


aE 

®xx  • ®yy.** 


^ [-»4s 


(89) 


% » 

4 1 

ahoira  that  alto  factors  out  of  each  of  these  terms. 

Hehctf  the  theraial  streaaea  are  teen  to  )>e  directly  propor- 


t I 
<1, 
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X 10*?  (DiMnsionl««s) 


7ig«  6 trmximt  9htm«l  Stress  Distributions  tor  One* 
Dimensibnsi  Axisl  Heat  Flat  (a024*t3  Aluninum) 


ticmal  to  tho  flux  density*  #|*  sad  the  plate  Sickness*  1. 

Mhen  expresses  in  s’oonsistent  set  of  units*  the  psraneter 
o 

h»9  the  inverse  dieSnsions  of  thermal  diffusivity. 
Bence*  nultiplying  both  sides  of  eqn  (89)  1^  results 

in  a ncm*dlmensionalised  stress  paraaeter 


Oxxic 


d£ 


/ T*ds^ 


^5E. 


h 


2h»  -h 


h 

/ T*sds 


] 


(90) 


where 


It  should  be  noted  that  eqn (90)  and  eqn  (91)  are  expressed 
in  terms  of  the  thermal  stress  and  heat  conduction  coordi- 
nate  systems*  respectively.  Eqn  (90)  is  plotted  in  Fig.  6 
for  2024»T3  aluminum.  If  the  material  property  parameters 
outside  the  brackets  in  eqn  (90)  were  transposed  to  \:he 
other  side  of  the  equatidn*  another  dimensionless  stress 
parameter  would  be  defined*  which  is  independent  of  material 
within  the  general  class  (hcxttogeneoUs*  isotropic)  of  interest. 

In  Pig.  6 the  curve  for  ict/l*  » 1.142  represents  the 
steady  state  stress  distribution  which  remains  constant  for 
all  rt/i^  > 1.  This  curve  is  the  same  as  that  shown  in  Pig. 

5 idiich  is  based  on  the  analytical  approximation  given  by 
eqn  (79).  For  earlier  tines*  the  develointent  of  the  thermal 


{ 
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•tresses  is  sIkmui  by  the  other  curves  in  rig.  6.  Initially, 
as  shown  by  the  curve  for  rt/i*  • .00571,  only  the  plate 
voluas  near  the  heated  surface  is  heated  appreciably  and  re« 
latively  large  coapressive  stresses  are  quickly  generated 
at  this  surface.  As  the  teaperature  drops  to  near  zero  be« 
yood  this  heated  layer,  the  esq^sion  of  the  heated  layer 
creates  tensile  stresses  in  the  unheated  area  idiich  reach 
a peak  near  the  front  of  the  relatively  unheated  zone,  then 
gradually  decay  becoaing  cosq>ressive  again  at  the  unheated 
surface. 

Both  the  cospressive  and  tensile  maxisnjzi  stresses  in- 
crease with  tine,  attaining  peak  values  prior  to  reaching 
the  steady  state  values.  The  peak  tensile  stress  is  only 
slightly  larger  than  the  steady  value,  which  occurs  at 
s/i  « 0.5.  However,  the  peak  cooq>ressive  stress  (at  the 
heated  surface)  is  approximately  20  percent  larger  than  the 
steady  state  value.  Fig.  6 also  shows  that  the  peak  com- 
pressive stress  occurs  at  an  earlier  time  than  the  peak 
tensile  stress. 

The  tenoral  development  of  the  peak  tensile  and  com- 
pressive stresses  is  shown  in  Figure  7.  The  vertical 
dashed  line  at  Kt/i*-  ■ 4.67  represents  the  time  at  which  the 
teaq>erature  variation  through-the-thickness  of  the  plate  be- 
comes less  than  1/10  that  of  the  heated  surface  twnperature. 
This  in  arbitrarily  taken  as  an  appriximation  of  the  time 
when  the  flux  vector  is  no  longer  primarily  axial  since  the 
radial  thermal  gradients  at  the  outer  diameter  of  the  beam 


would  bw  wsq^ectwd  ijo  be  large  relative  to  the  exiel  grad- 
imits  tmder  the  bemm. 

The  Blight  jr^eak  in  the  tensile  stress  near  Kt/t*  • 

0.1  is  iq»parent.  The  peaking  of  the  cosg>ressive  stress  is 

seen  to  be  wuch  wore  significant  and  extending  over  a larger 

period' of  tine.  It  is  concluded  that  the  naxianm  themal 

stresses  producedl  by  one'dinensional  axial  heating  occur 
« 

just  prior  to  attaining  steady  state  heat  conduction,  and 
that  these  Stresses  are  defined  by  the  following  equations 
as  derived  fron'Fig.  7. 


(Oxx<  \ 

W.t  / 


* 3.30  X 10“*0ict/t*  » o.r 


/ ^xx*^  \ 

W.t  / 


max 

tensile 


» -.7.04  x 10*'*0ict/l*  - 0.06 

max. 

coB^ressive 


To  assess  the  relative  importance  of  the  melting  and 

•5 

iherxial  degradation  damage  mechanisms  for  comparison  with 
thermal  stresses*  for  the  axial  flux  ckse«  the  times  required 
for  complete  melt**thro;^gh^  and  for  cgmplete  thermal  degrada- 
tion ai'e  useful.  The  time  required  for  complete  melt- 
through,  assuning  one>dimei^lonal  astial  flux,  temperature 
independent  properties,  and  immediate,  melt  removal,  , is 
given  by  (kef,.  19 $2) 


t,  - |i-tCp(T«-»ol+I»l 


fiM  defined  in  ^ectiqn  VX,  thermal  degradation  le  a 


o 


fimeti<»i  of  tMporatore  only*  with  •sgentlally  zero  strength 
resMiinifig  if  local  tenperatores  exceed  370*C.  Analogous  to 
the  one*^iaensional  axial  flux  ■alt«>through  tiae»  tt,  a 
heating  tins  for  coeplete  themal  degradation  can  be  defined 

by 

to  - {i-CpCTo-To)  (95) 

idiera 

to  * heating  tine  for  eong>lete  themal  degradation 
of  structural  strength  (sec) 

Tq  • tesg>eratare  above  which  material  has  essen~ 
tially  zero  strength 
- 370»C 


Dividing  by  t|  gives 

*2  Tq-Tq 
t,  * Ti5=5p+iVCp 

It  is  seen  that  to  is  always  less  than  ti  for  structural 
Bietals  since  is  generally  greater  that  Tq.  For  2024-T3 
aluminum,  assuming  « 0, 


to  **  0.43  ti 


(97) 


Herce,  complete  thermal  degradation  of  the  structural 
strength  properties  of  the  plate  would  be  expected  to  occur 
prior  to  ccmiplete  melt-through. 

Although  not  directly  pertainent  to  the  present  study, 
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tlMi  bMtlng  tiMi  to  coapleta  thoraal  degradation  for  other 
ctriictaral  aateriale  wae  investigated  briefly.  Based  on 
extrapolation  of  the  tbexaal  degradation  data  in  Ref  9,  the 
folloifing  table  presents  t^/ti  as  calculated  from  eqn  (96). 
fhersMil  properties  are  talcen  froai  Ref  32. 


Material 

T«(«C) 

*PfC) 

gp.U/ffl.*C) 

AS-31B 

Nagnesiim 

605 

425 

338 

1.21 

0.48 

301 

Stainless 

Steel 

1400 

870 

290 

0.42 

0.42 

6At-4V 

Titanium 

1600 

870 

390 

0.77 

0.41 

c 


Bence,  it  appears  that  c<»iplete  thermal  degradation  precedes 
melting  for  most  structural  metals. 

A qualitative  comparison  of  the  thermal  degradation  and 
thermal  stress  damage  mechanisms  for  the  axial  flux  case  can 
be  obtained  by  examining  the  ratio  of  heating  time  to  com- 
plete thermal  degradation  tim«.  and  comparing  it  to  the  ratio 
of  maximum  compressive  thermal  stress  to  room  temperature 
conqpressive  yield  strength.  This  oonqmrison  is  made  in  Fig. 
8 for  a particular  value  of  absorbed  heat  flux  density. 

Since  both  of  these  dimensionless  ratios  are  linear  func- 
tions of  the  absorbed  heat  flux  density  (see  oqns  (90) 
and  (95)),  changes  in  this  parameter  do  not  affect  the  re- 
lative positions  of  the  tvo  curves. 

The  solid  curves  ai^ly  v^n  xt/i*  * 'O.OS  which  corres- 
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—mKt/t^  m 0.06 
• 0.01 


Flj$.  8 Comparisons  of  «Peak  Thermal  Stresses 
and  Thermal  Defp^adation  Heating  Time 
for  One^Dimensional  Axial  Heat  Flux 
C2024-T3  Aluminum) 
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•%»'  •* 


ponds  to  tho  MxiJMMi  value  of  the  peak  ox^ressive  stress 
(see  Fig.  7).  The  stress  ratio  is  cgsputed  fron  eqn  (93) 
by  substituting  the  Material  property  values  for  2024~t3 
aluMlnuM  and  » 1500  v/cm*  which  is  approximately  tho 
MSxiMUM  value  of  interest  for  this  study.  The  fo^ldi^ng 
relationship  is  obtained. 


/ff—  \ 

l?co  j, 


* 0.9e9(l/cB)i0xtA*  * 0.06 


(98) 


Similarly,  eqn  (9!  gives 


1.39  (ccK/sec)^ 


(99) 


,,---ncroM  Kt/t*  >•  0.06, 


t » O.llTlMsec) 


(100) 


4 

Substituting  (100)  into  (99)  results  in 

« 0.162 (l/cm)i  (101) 

ta 

Bgns  (98)  a*;:!  (101)  are  plotted  as  solid  curves  in  Fig.  8. 
The  initial  conclusion  drawn  from  this  comparison  is  that 
the  damage  fron  these  two  sources  is  of  the  same  order  of 
magnitude  at  times  defined  by  Kt/t*  ■ 0.06.  This  is  a quail 
tatlve  conclusion  in  that  peak  thermal  stresses  represent 
the  onset  of  material  failure  rather  than  through- the> thick- 


ness  failure  represented  by  the  thermal  degradation  ratio. 
Even  though  the  stress  ratio  exceeds  the  thc<.mal  degrada- 
tion ratio  by  about  a factor  of  six,  consideration  of  the 
additional  time  required  to  completely  fail  the  plate 
through-the-thickness  by  both  compressive  and  tensile  ther- 
mal * stresses  would  likely  compensate  for  this  difference 
and  might  reverse  the  relationships. 

For  comparison  to  the  maximum  compressive  stress,  the 
corresponding  stress  ratio  for  maximum  tensile  thermal 
stresses  is  shown  on  Fig.  8 as  derived  from  eqn  (92).  Ten- 
.«ile  stresses  are  concluded  to  be  a relatively  insignificant 
contributor  to  damage  for  axial  flux  heating. 

Also  shovm'in  Fig.  8 are  the  corresponding  relation- 
ships for  an  earlier  time,  <t/l^  « 0.01.  In  this  case  the 
stress  ratio  is  determined  from  Pig.  7.  Since  t/tj)  is 
linear  in  time  and  the  peak  compressive  stress  is  a loga- 
rithmic function  of  time,  the  stress  ratio  exceeds  the  ther- 
mal degradation  ratio  by  an  even  larger  factor.  Consequen- 
tly the  relative  significance  of  thermal  degradation  in- 
creases with  time.  Since  the  heating  conditions  for  the 
solid  curves  in  Fig.  8 generate  the  maximum  thermal  stresses 
for  the  range  of  beam/target  parameters  considered,  it  is 
seen  that  the  maximum  value  for  the  compressive  stress  ratio 
is  about  0.65  (at  the  maximum  plate  thickness  of  interest, 
t XR  0.65  cm) . 

Hence,  for  the  axial  flux  conditions  considered,  peak 
thermal  stresses  do  not  reach  material  failure  values. 
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Connie l:c  thermal  degradation  occurs  at  t/£  « 0.719  sec/cm 
as  determined  from  eqn  (99) . This  corresponds  to 


Kt 


0.370 


(102) 


^ich  varies  between  0.57  and  3.70  for  the  range  of  plate 
thicknesses  considered  (0.1>l>0.65cm) . That  is,  complete 
thermal  degradation  does  occur  for  the  axial  heat  flux  case 
within  the  range  of  beam/ target  parameters  of  interest. 

Since  the  range  of  Kt/l^  in  which  complete  degradation  can 
occur  exceeds  the  ict/i^  values  for  maximum  thermal  stresses, 
the  actual  failure  mechanism  is  expected  to  be  the  combined 
actions  of  thermal  stresses  and  thermal  degradation.  How> 
ever,  since  the  majority  of  the  plate  thickness  under  the 
beam  is  subjected  to  tensile  stresses  (Pig.  5),  and  these 
tensile  stresses  do  not  become  significant  relative  to  room 
temperature  ultimate  tensile  strength  for  the  range  of  heat- 
ing parameters  considered,  thermal  degradation  will  contri- 
bute much  more  to  complete  through-the-thicknoss  failure 
than  thermal  stresses.  It  is  concluded  that  the  heating 
time  to  complete  thermal  degradation,  tjj,  provides  a rea- 
sonable approximation  to  the  time  required  for  through- the- 
thickness  material  failure  under  the  combined  actions  of 
thermal  stresses  and  degradation  for  the  case  of  axial  flux 
heating. 

Comparisons  of  thermal  stresses  generated  during 
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one>dimen8ional  axial  heat  conduction  to  those  associated 
with  one**diinensional  radial  conduction  are  discussed  in 
later  paragraphs  of  this  section. 

Radial  Flux  Analyses 

The  numerical  radial  flux  computer  program  is  based  on 
the  numerical  thermal  model  described  in  Section  IV.  The 
total  progr^un  is  nzuned  TSTRESS  and  is  listed  in  Appendix  A 
along  with  sample  program  outputs.  The  radial  thermal  model 
described  in  Section  IV  is  subroutine  TEMPRAD  in  Appendix  A. 
Based  on  the  radial  temperature  distributions  output  from 
TEMPRAD,  numerical  integration  of  thermal  stress  equations 
(55)  & (56)  is  accomplished  by  subroutine  STRESSR  to  deter- 
mine the  radial  flux  stress  components.  A separate  algo- 
rithm (see  Appendix  B)  was  developed  to  compute  stresses 
near  r «*  0 as  equations  (55)  & (56)  are  indeterminate  there. 

The  one-dimensional  radial  heat  conduction  model  im- 
plies that  the  plate  is  circular  in  planform.  However,  for 
those  cases  where  the  plate  diameter  is  large  enough  that 
the  plate  size  no  longer  significantly  affects  the  heat 
conduction  or  thermal  stresses,  the  plate  can  be  considered 
to  be  rectangular  with  width  equal  to  the  circular  plate 
diameter.  The  length  of  the  rectangular  plate  can  be  con- 
sidered to  be  any  length  equal  to  or  greater  than  the  cir- 
cular plate  diameter.  The  residual  tensile  strength  of  the 
rectangular  plate  in  the  lengthwise  direction  can  then  be 
used  as  a measure  of  the  damage  inflicted  by  the  mechanisms 


of  thernal  degradation  of  structural  properties,  themal 
stresses,  and  nelting. 

Subroutine  STRONG  conputes  the  radially  averaged  re» 
sidual  tensile  strength  of  the  plate  (RESDEG)  due  tc  ther- 
mal degradation,  based  on  the  radial  temperature  distribu- 
tion and  the  half-hour  exposure,  thermal  degradation  of  ul- 
timate tensile  strength  data  for  2024-T3  aluminum  given  in 
Ref  9.  These  data  are  approximated  by  the  analytical  func- 
tions given  in  Section  VI.  This  subroutine  provides  a mea- 
sure of  the  damage  caused  by  thermal  degradation  of  tensile 
strength. 

Other  damage  measures  are  computed  in  subroutine 
RESTRNG,  Residual  tensile  strength  due  to  thermal  stresses 
(RESTRSS)  is  computed  based  on  thermal  tensile  stress  dis- 
tribution by  assuming  ultimate  tensile  stress  remains  at 
the  room  temperature  value.  Residual  strength  due  to  the 
combined  action  of  thermal  stresses  and  thermal  degradation 
is  computed  as  RESTRS2.  Structurally  failed  areas  due  to 
thermal  stresses  alone  (PPS),  thermal  degradation  alone 
(RFD) , and  in  combination  (RFDS)  are  calculated  in  subrou- 
tine RADFAIL.  Failed  radius  due  to  melting  a hole  (RMELT) 
is  computed  in  the  main  program  (beginning  at  line  153) . 

A parallel  calculation  of  the  effect  of  stress  concen- 
tration is  made  in  subroutine  RAOSC.  This  calculation  is 
based  on  holes  caused  by  thermal  stresses  and  degradation 
in  combination  (RFDS)  as  it  is  assumed  that  hole  stress 

f 

concentrations  associated  with  melted  holes  are  negligible 
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due  to  relief  provided  by  plastic  flow  associated  with 
C ^ high  temperatures. 

Structural  failure  may  be  caused  by  thermal  stresses 
or  degradation,  or  their  combined  actions.  Two  alternative 
treatments  of  structurally  failed  areas  were  used.  In  one 
case,  the  failed  area  is  assumed  to  be  immediately  removed 
fr<»n  the  plate,  thus  forming  a hole.  In  the  other  case, 
the  failed  area  is  assumed  to  remain  in  the  plate  and  con- 
tinue to  absorb  and  conduct  heat. 

In  the  case  where  the  failed  area  is  retained  in  the 
plate,  melting  is  the  only  mechanism  which  can  generate  a 
hole.  At  each  time  increment,  theirmal  stresses,  residual 
strengths  and  failed  areas  are  all  computed  based  on  the 
f assumption  that  the  owrrent  failed  area  has  zero  strength, 

V 

that  is  only  the  temperature  distribution  outside  the  failed 
area  is  considered.  The  temperature  distribution  is  based 
on  the  heat  flux  into  the  plate  outside  the  melted  hole. 

The  failed  area  V7ill  either  be  equal  to  or  greater  than  the 
melted  hole. 

In  the  case  where  the  failed  area  is  assumed  to  be  im- 
mediately removed  (by  an  unspecified  source  such  as  aero- 
dynamic pressure) , the  actual  hole  size  at  any  time  is  the 
larger  of  the  melted  hole  and  the  failed  area.  Temperatures, 
thermal  stresses,  residual  strengths  and  failed  areas  are 
computed  based  on  this  actual  hole  size. 

Fig.  9 (solid  curves)  shows  a typical  comparison  of 
C)  the  respective  hole  sizes  for  the  above  described  ilterna- 
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tive  failure  nodes.  The  plots  are  based  on  data  output 
from  the  TSTRBSS  program  at  one  second  time  intervals. 

The  data  are  connected  by  straight  lines  which  results  in 
the  incr^oental  characteristic  of  the  plots.  The  bottom 
solid  curve  is  for  the  case  in  which  only  melting  is  consid- 
ered to  form  a hole  in  the  plate.  For  these  conditions  a 
failed  area  is  generated  by  the  combined  actions  of  thermal 
stresses  and  thermal  degradation  and  the  radius  of  this 
failed  area  is  given  by  the  dashed  curve.  The  failed  area 
is  assumed  to  remain  in  position  in  the  plate.  For  the  al- 
ternative situation  in  which  the  failed  area  is  assumed  to 
be  immediately  removed  from  the  plate#  the  resulting  hole 
radius  is  given  by  the  upper  solid  curve  in  Fig.  9.  These 
results  are  typical  for  the  range  of  parameters  studied  in 
that  the  structurally  failed  area  due  to  the  combined  action 
of  thermal  stresses  and  degradation  is  always  larger  than 
holes  formed  by  melting#  and  that  the  failed  area  is  gener- 
ated prior  to  melting  a hole. 

It  should  be  noted  that  although  the  one-dimensional 
radial  flux  code  is  strictly  applicable  only  when  the  dimen- 
sionless power  per  unit  thickness#  Pj^^#  is  less  than  5#  much 
of  the  data  presented  is  based  on  heating  conditions  for 
which  > 5.  This  is  done  to  facilitate  data  presentation 
as  the  trends  shown  are  the  same  for  both  conditions,  but 
much  longer  computer  run  times  are  required  to  obtain  struc- 
tural failure  and  melting  effects  when  P^^^  < 5. 

Fig.  10  is  a typical  plot  of  principal  thermal  stress 
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RADIUS  (CM) 

Fig.  10  Thermal  Stress  Distributions.  1-D  Radial  Heat  Flux,  Gaussian  Heat  Flux  Density 


distributions  for  a solid  disk  at  some  time  prior  to  melting 
or  structural  failure.  An  understanding  of  thermal  stress 
distributions  prior  and  subsequent  to  hole  formation  is 
necessary  to  explain  the  relative  damage  modes.  Fig.  10 
shows  that  the  radial  (a,.)  and  tangential  (Cg)  stresses  are 
compressive  near  the  plate  and  beam  center.  The  radial 
stress  becomes  less  compressive  with  increasing  radius, 
becemting  zero  at  the  plate  edge.  The  tangential  stress  be- 
ccxnes  tensile  just  outside  the  heated  cylinder,  rapidly 
reaching  a peak  value  before  gradually  decreasing  to  zero 
at  the  plate  edge. 

The  maximum  shear  stress  is  equal  to  the  one-half  the 
maximum  difference  between  any  two  of  the  three  principal 
orthogonal  stresses.  In  the  case  of  plane  stress,  the 
stress  component  normal  to  the  plane  of  the  plate  is  zero. 
Inside  the  heated  cylinder  the  maximum  difference  in  stress 
components  is  between  the  radial  and  normal  componetjts. 
Outside  the  heated  core,  the  maximum  shear  stress  is  deter- 
mined by  the  difference  between  the  radial  and  tangential 
components.  This  change  in  the  maximum  shear  stress  deter- 
mination accounts  for  the  inflection  In  the  maximum  shear 
stress  distribution  shown  in  the  figure. 

A typical  temperature  distribution  at  some  time  prior 
to  formation  of  a hole  is  illustrated  in  Fig  11.  Both  melt- 
ing and  thermal  degradation  are  direct  functions  of  tempera- 
ture (see  Section  VI),  and  knowledge  of  the  general  form  of 
the  temperature  distribution  aids  In  understanding  the  reXa- 
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tive  influence  of  these  damage  mechanisms.  The  figure  shows 
the  typically  rapid  decrease  in  temperature  within  the  beam 
radius  for  a Gaussian  distribution  heat  flux  density. 

Pig.  12  presents  the  residual  tensile  strength,  nor' 
malized  by  the  unheated  ultimate  tensile  strength,  as  a 
function  of  heating  time  for  each  of  the  damage  mechanisms. 
In  this  case,  the  failed  area  is  assumed  to  be  immediately 
removed,  thus  forming  a hole.  Since  the  failed  area  radius 
is  always  equal  to  or  greater  than  any  melted  hole  radius, 
no  damage  due  to  melting  occurs. 

The  effect  of  assuming  that  the  failed  area  remains  in 
the  plate  but  does  not  transmit  stresses  and  that  only  melt- 
ing generates  a hole  is  shown  in  Fig.  13.  Comparison  with 
Pig.  12  shows  that  the  assumption  that  the  failed  area  re- 
mains in  place  results  in  slightly  greater  strength  reduc- 
tion, again  due  to  the  larger  heat  absorbing  surface  for 
this  situation.  Both  figures  show  that  the  largest  contri- 
butor to  tensile  strength  reduction  is  the  structurally 
failed  area  arc"  that  thenial  stresses  contribute  more  than 
thermal  degradation.  The  respective  radii  of  the  melted 
hole  and  failed  area  for  the  heating  conditions  of  Fig.  13 
are  given  by  the  solid  curves  in  Pig.  14.  The  dashed  curve 
is  for  the  alternative  assumption  that  the  failed  area  forms 
a hole. 

The  effect  of  varying  laser  beam  diameter  on  residual 
tensile  strength  is  shown  in  Fig.  15  for  the  situation  of 
retained  failed  areas.  For  smaller  beam  diameters  the  peak 
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heat  flux  density  is  larger  resulting  in  earlier  structural 
failure.  As  the  failed  area  diameter  approaches  the  beam 
diameter,  the  rate  of  strength  reduction  becomes  increasing- 
ly slower  because  the  failed  area  is  the  primary  factor  in 
determining  residual  tensile  strength.  Residual  strength 
continues  to  decrease  until  the  melted  hole  approaches  the 
beam  size  after  which  heat  absorption  by  the  plate  ceases. 

The  larger  beam  diameter  effectively  postpones  struc- 
tural failure  and  melting  until  a larger  area  of  the  plate 
is  heated  significantly  resulting  in  strength  reduction  over 
a larger  portion  of  the  plate.  As  shown  in  Fig.  16,  the 
failed  area  grows  more  rapidly  for  the  larger  beam  diameters 
eventually  becoming  larger  than  for  the  smaller  beam  dia- 
meters. 

The  effect  of  beam  power  on  residual  strength  is  shown 
in  Fig.  17.  The  beam  powers  selected  for  this  comparison 
give  peak  absorbed  flux  densities  corresponding  to  those 
for  the  comparison  of  beam  diameter  effects  in  Fig.  15. 

The  corresponding  failed  area  radii  are  compared  in  Fig.  18. 
Increasing  beam  power  is  seen  to  decrease  the  time  required 
to  achieve  structural  failure  and  melting,  and  to  generate 
larger  failed  areas  (and  melted  holes)  and  reduce  residual 
tensile  strength  at  any  given  time. 

Fig.  19  shows  the  effect  of  plate  thickness  on  residual 
strength.  The  effects  are  similar  to  those  for  beam  power 
with  more  rapid  strength  reduction  consistently  occur ing 
for  ♦•h inner  plates. 
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17  Effect  of  Atsorbed  Power  on  Residual  Tensile  Strength 
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Pig.  18  Effect  of  Absorbed  ?ov/er  on  Melted  Hole  and  Structurally 
Failed  Area  Radii 


Pig.  19  Effect  of  Target  Plate  Thickness  on  Residual  Tensile  Strength 


One  objective  of  the  present  study  was  to  determine  the 
relative  importance  of  the  thermal  stress  damage  mechanisiit 
relative  to  melting  and  thermal  degradation.  From  the  data 
presented  thus  far  for  the  case  of  ortc>diroensional  radial 
fluxr  one  major  conclusion  is  that  the  primary  mode  of  dam- 
age (measured  in  terms  of  residual  strength)  is  the  genera- 
tion of  a structurally  failed  area  by  the  combined  action 
of  thermal  stresses  and  degradation.  The  relative  contribu- 
tions of  each  of  these  mechanisms  to  structural  failure  is 
thus  of  interest.  A qualitative  assessment  of  the  relative 
contributions  can  be  obtained  by  looking  at  the  peak  stress- 
es and  temperatures  (outside  the  failed  area)  as  shown  in 
Fig.  20.  This  plot  is  typical  for  all  the  cases  investigat- 
ed in  that  once  a structur  lly  failed  area  is  initiated,  the 
tangential  compressive  stresses  remain  at  relatively  high 
levels  while  the  other  stress  components  are  reduced  signifi- 
cantly. Fig.  21  and  Fig.  22  show  that  peak  tensile  stresses 
remain  very  small  for  a wide  range  of  beam  powers  and  dia- 
meters. Since  both  peak  tangential  compressive  stress  and 
maximum  temperature  (outside  failed  area)  occur  at  the  radius 
of  the  failed  area,  it  is  concluded  that  these  are  the  mech- 
anisms v/hich  continue  to  expand  the  failed  area  with  increas- 
ing time.  Tl'at  is,  the  failed  area  expands  due  to  compres- 
sive yield  failure. 

Although  thermal  degradation  is  not  generally  a linear 
function  of  temperature,  for  the  analytical  approximations 
used  in  this  study  (Section  VI)  and  the  range  of  temperatures 
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outside  the  failed  area  shown  in  Fig.  20,  the  relationship 
is  essentially  line2ur.  Consequently,  the  temperature 
curve  in  Pig.  20  is  a fairly  good  approximation  of  the  com- 
pressive strength  reduction  due  to  thermal  degradation. 

The  general  conclusion  is  then  made  that  the  compressive 
tangential  thermal  stresses  contribute  significantly  more 
to  structural  failure  than  does  thermal  degradation. 

A quantitative  assessment  of  the  relative  contributions 
to  structural  failure  does  not  seem  feasible.  The  essential 
conclusion  is  that  qualitatively,  compressive  tangential 
thermal  stress  at  the  edge  of  a failed  area  is  the  primary 
contributor,  along  with  theimal  degradation,  to  the  expan- 
sion of  the  deimaged  or  failed  area. 

Prom  the  analysis  of  the  thermal  stresses  for  the  one- 
dimensional radial  flux  code  (reference  sample  output  in 
Appendix  A),  it  was  found  that  prior  to  structural  faiJurc 
or  melting,  thermal  stress  is  a linear  function  of  the  beam 
power  and  consequently  linear  in  absorbed  power,  peak  ab- 
sorbed flux  density,  and  peak  absorbed  fluence  through  the 
following  relationships. 


Pa  = 

(103) 

Pa  - 

(104) 

Ppa  • 

^pa^ 

(105) 

where 
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Thermal  stresses  also  were  found  to  be  inversely  proportion- 
al to  the  target  plate  thickness.  These  empirical  findings 
were  verified  by  examination  of  the  algorithms  used  for 
temperature/  eqn  (27),  and  thermal  stresses,  egns  (55)  and 
(56) . These  relationships  suggested  two  dimensionless 
stress  parameters  which  were  studied  to  see  whether  they  had 
any  particular  invariant  qualities  which  might  serve  to 
characterize  the  thermal  stress  problem. 

The  parameter  ^ , where  of  represents  any  thermal 

stress  component,  expressed  in  consistent  units,  has  the 
units  of  sec/cm*  which  are  the  inverse  units  of  thermal 
diffusivity,  x.  Hence,  the  parameter  is  dimension- 

less. The  parameter  2^  has  the  units  of  seconds. 

ipa 

Dividing  by  the  heating  time  is  equivalent  to  replacing  the 
peak  absorbed  flux  density  in  the  denominator  by  the  peak 
absorbed  fluence,  giving  another  dimensionless  parameter 
Fig.  23  compares  this  parameter,  using  peak  ten- 
sile stress,  for  three  beam  diameters,  as  a function  of 
heating  time. 

The  peak  absorbed  fluence  used  in  Fig.  23  is  that  for 
a continuously  solid  plate  without  regard  to  the  effect  of 


al 

^a 
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the  expanding  hole  in  the  plate,  ffhen  this  effect  is  ac- 
counted for  by  uuing  the  actual  peak  absorbed  fluence,  the 
results  are  less  orderly  than  in  Fig.  23  due  to  the  step 
function  effect  of  the  hole  growth  in  the  numerical  model. 

Fig.  24  shows  the  other  dimensionless  parameter  des- 
cribed above  for  the  same  conditions  of  Fig.  23.  Compari- 
son of  the  two  plots  shows  that  there  is  less  spread  in  the 
data  based  on  , and  that  the  relative  positions  of 

the  curves  based  on  is  more  indicative  of  the  rela- 

?ir 

tive  peak  tensile  stresses,  being  a linear  function  of  the 
peak  tens^e  stress  (see  Fig.  21) . 

Figure  25  shows  the  development  of  peak  thermal  stresses 
for  a typical  case  in  v;hich  the  heat  conduction  is  initially 
axial.  The  figure  shows  the  theoretical  stresses  for  both 
one-dimensional  axial  and  radial  conductivity  models  for 
identical  conditions.  Radial  flux  stresses  are  given  for 
a uniform  beam  flux  disitribution.  The  uniforja  beam  is  con- 
sistent with  the  constant  flux  distribution  assumed  in  the 
axial  conduction  model.  Fig.  25  shows  that  the  initial 
transient  and  subsequent  steady  state  peak  stresses  for 
one-dimensional  axial  flux  remain  fairly  small  relative 
T.o  compressive  yield  ot  ultimate  tensile  strengths  (23,  ) 

n/cm*  and  44,820  n/cm* , respectively). 

For  the  flux  density  and  thickness  listed  on  Fig.  25, 
the  initial  conduction  is  primarily  axial  through-the- 
thickness  if  the  beam  diameter  is  at  least  5.87  cm.  In 


this  case,  Pjj,a  ^ 70.0  as  determined  from  eqn  (75).  Hence, 
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heat  flux  for  this  case  would  remain  axial  until  a hole  is 
melted  through  the  plate , No  appreciable  radial  flux  would 
develop  in  this  time.  Consequently,  the  radial  flux  results 
given  in  Fig.  25  are  strictly  academic  and  serve  here  to 
indicate  the  trends  in  thermal  stress  development  for  those 
cases  where  the  flux  vector  becomes  predominatly  radial  be- 
fore melting  occurs. 

These  results  illustrate  a bound  on  the  radial  flux 
model.  If  ^ 70  at  the  onset  of  heating,  the  target 
plate  will  melt  through  without  any  transition  to  a radial 
flux  case.  Hence,  the  radial  flux  model  is  not  generally 
applicable  if,  initially,  noted  earlier  in 

this  section,  if  P£-^  < 5 at  the  onset  of  heating,  the  heat 
flux  vector  is  primarily  radial  and  remains  so  as  long  as 
heating  continues.  In  Intermediate  cases  (5<P;,a<70) , the 
flux  vector  would  be  expected  to  become  increasingly  radial 
at  long  enough  heating  times.  From  eqn  (87),  the  flux  vector 
would  be  expected  to  be  primarily  radial  at  lea.-3t  by  the 
time  Kt/i*  is  about  5.  The  flux  vector  becomes  primarily 
radial  at  ?yer  shorter  heating  times  for  decreasing  P^a* 

The  initial  flux  vector  would  be  expected  to  be  predominant- 
ly radial  for  values  of  p£jj  near  5 and  to  be  mostly  axial 
for  values  of  approaching  70.  The  axial  flux  thermal 
stress  model  is  appli<iable  for  all  heating  times  through 
melting  if  p£j^  > 70. 

Correlation  of  dimensionless  power  per  unit  thickness, 
p£j^,  with  the  time  to  initial  structural  failure  as  predict- 

Q6 


•d  by  the  one^dimensional  radial  flux  model  was  investigated 
briefly.  Fig.  26  presents  the  results  of  this  correlation, 
using  data,  denoted  by  the  circluar  data  points,  that  had 
been  developed  previously  as  part  of  the  other  analyses  con- 
ducted. 

A definite  correlation  is  apparent  which  also  appears  ) 
to  be  a function  of  the  laser  beam  diameter.  The  data 
presented  include  all  of  the  pre-existing  data  available 
trhich  was  applicable.  That  is,  no  sorting  of  the  data  was 
done  to  bias  the  results.  As  noted  previously,  most  of  the 
data  frexn  the  radial  flux  code  are  for  heating  conditions 
with  P£a  > 5 which  violates  the  assumed  bound  on  the  appli- 
cability of  the  model.  Consequently,  additional  analyses 
are  required  to  extend  the  correlations  to  values  of  < 

5.  The  validity  of  the  radial  code  correlations  for  P|^a  ^ 

5 were  qualitatively  assessed  by  plotting  corresponding 
failure  times  for  ona-dimensional  axial  flux  > 70)  on 

Fig.  26. 

For  one-dimensional  axial  flux,  the  heating  time  to 
complete  thermal  degradation  can  be  expressed  as  follows  by 
combining  eqns  (75)  and  (95)  and  substituting  tns  2024-T3 
material  values. 


tf)  «>  0.658  (D  in  cm)  (106) 

This  equation  is  plotted  in  Fig.  26  for  three  diameters  cor- 
responding to  the  radial  flux  data.  The  radial  and  axial 
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flux  data  are  not  directly  dcxnparable  because  the  radial 
data  are  based  on  a Gaussian  heat  flux  density  distribution 
and  the  axial  model  assximes  a uniform  distribution.  To 
make  a more  consistent  comparison,  four  radial  flux  code 
cases  with  a be£un  diameter  of  10  cm  were  analyzed  for  a uni- 
form flux  distribution.  These  results  are  denoted  in  Fig. 

26  by  the  x-symbol.  The  results  of  this  limited  comparison 
indicate  that  a possibly  consistent  relationship  between 
failure  time  and  may  exist  over  a wide  range  of 
values.  Further  study  of  this  relationship  is  necessary  to 
obtain  more  definite  conclusions. 

Correlations  of  with  the  size  of  structurally  fail- 
ed damage  areas  as  a function  of  the  laser  beam  diameter  are 
shown  in  Fig.  27  for  the  cases  listed  in  the  previous  fig- 
ure. The  da\;a  plotted  are  the  radii  of  the  structurally 
failed  area  at  a heating  time  of  15  seconds.  Smooth  curves 
were  drawn  through  the  data  points. 

It  was  concluded  that  Pj^^^  does  provide  a consistent 
correlation  with  the  structural  damage  size.  Quantitative 
definition  of  this  correlation  was  not  established  and  is  a 
topic  requiring  additional  study.  The  correlation  examined 
here  is  actually  with  the  ”dimensional"  power  per  unit  thick- 
ness, since  the  additional  quantities  in  P^^  are  mat- 

erial properties  and  only  a single  material  was  considered. 
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VIII.  Conclusions  and  Recommendations 


Conclusions 

For  both  limiting  cases  of  one-dimensional  axial  cind 
radial  flux,  the  primary  mode  of  damage  is  by  generating 
a failed  area  under  the  laser  beam.  For  axial  flux  condi- 
tions, the  primary  damage  mechanism  is  thermal  degradation 
of  stru'ctural  strength,  based  on  the  half-hour  heating  ex- 
posure data  of  (Ref  9) . For  the  range  of  beam/target  para- 
meters considered  here  (see  Section  II),  thermal  stresses 
are  relatively  small.  For  the  case  of  one-jimensional 
radial  flux,  the  primary  damage  mode  is  structural  failure 
due  to  the  combined  actions  of  thermal  stresses  and  thermal 
degradation.  If  the  structurally  failed  area  is  retained 
in  the  plate,  melting  may  occur,  but  never  precedes  struc- 
tural failure.  Detailed  conclusions  are  given  in  the  fol- 
lowing paragraphs. 

• Axial  Flux.  For  this  study,  axial  flux  was  defined  to 
exist  if  the  dimensionless  power  per  unit  thickness,  Pji^r 
is  greater  than  about  70.  The  flux  was  defined  to  become 
primarily  radial  for  dimensionless  thermal  times,  tct/A®, 
greater  than  about  5 since  at  this  time  the  through— the- 
thickness  temperature  variation  becomes  less  than  ten  per- 
cent. TWO  time  regimes  Were  identified  as  significant  in 
characterising  the  thermal  stres$es  for  axial  flux.  The 
maximum  thermal  stresses  for  the  beam/target  parameters  of 


interest  occur  during  a transient  stage  when  Kt/l^  < 1. 
These  maximum  stresses  become  constants  when  expressed  in 
the  following  dimensionless  form 


* 3.30  X 10“*§ict/i*  = 0.1  (107) 

/peak 

tensile 

( peak  » - 7.94  x lO^^eKt/t*  * 0.06  (108) 

\ ^ /compressive 

The  peak  through-the-thickness  stresses  decrease  with  in- 
creasing time  until  they  reach  a steady  state  value  for  all 
Kt/l^  > 1.  For  these  steady  state  conditions,  the  peak 
thermal  stresses  are  slightly  less  than  the  maximum  values 
given  in  eqns  (107)  and  (108) . For  the  range  of  beam/tar- 
get parameters  considered  in  this  study,  the  thermal  stress- 
es remain  very  small  relative  to  target  material  strength. 

Since  both  melting  and  thermal  degradation  are  direct 
functions  of  temperature,  thermal  degradation  is  the  ini- 
tial damage  mechanism  because  the  temperature  necessary  to 
reduce  the  material  strength  to  zero  is  lower  than  the 
temperature  required  for  melting.  If  it  is  assumed  that 
any  material  failed  by  thermal  degradation  remains  in  place 
in  the  plate,  melting  will  subsequently  occur  with  contin- 
ued heating.  Any  material  melted  was  assumed  throughout 
this  study  to  be  immediately  removed. 

tot  axial  flux  conditions  the  damage  is  essentially 
restricted  to  the  area  under  the  incident  beam  spot.  Con- 
sequently, the  damage  size  vbuid  correspond  ctlosely  to  tne 

9i' 


team  diameter,  for  a spatially  uniform  flux  density  distri- 
bution. For  a Gaussian  beam,  the  damage  area  would  be  re- 
stricted to  an  area  near  the  team  axis  because  of  the  much 
lower  flux  density  near  the  beam  edges  which  would  tend  to 
develop  radial  flux.  The  time  required  to  fail  an  area 
equal  to  the  beam  spot  size  is  much  less  than  for  one-dimen- 
sional radial  flux. 


Radial  Flux.  The  heat  conduction  in  the  heated  target 
plate  was  defined  to  be  essentially  one-dimensional  radial 
if  < 5.  For  cases  where  5 < <70,  the  flux  would 

be  expected  to  become  increasingly  radial  with  increasirig 
heating  time.  The  nximerical  radial  flux  code  developed  as 
part  of  this -.study  would  be  expected  to  give  approximate 
results  for  heating  conditions  with  these  intermediate 
valued  of  Pj^  at  long  enough  heating  times  (Kt/Jl*>5)  and 
accurate  results-  for  p£J^  < 5.  Temperatures  and  thermal 
stresses  for  the  radial  code  are  compared  to  those  from  a 
two-dimensional  NASTRAN  finite  element  analysis  in  Appen- 
dix C. 


Based  on  the  data  from  the  6ne-dimensional  radial  anal- 

ysis,  it  was  concluded  that  the  primary  damage  mechanism 

for  radial  flux  heating  conditions  is  structural  failure  by 

, ^ 

thermal  stresses  and  strength  degradation  in  combination. 

The  primary  contributor  to  residual  tehsile  strength  reduc- 
tioh  is  from  the  reduction  in  net  section  load  carrying 
ability  resulting  frqm  &itructural  failure  of  a circular 
portion  of  the  plate.'  $tructural  failure  is  due  primarily 


! 


to  compressive  tangential  thermal  stresses » in  combination 
with  thermal  degradation.  Reductions  in  residual  strength 
due  to  thermal  degradation  and  thermal  stresses  outside  the 
failed  area  are  much  less  than  that  due  to  the  failed  area. 

Thermal  degradation  contributes  less  to  residual 
strength  reduction  outside  the  failed  area  than  do  thermal 
stresses  because  significant  temperatures  are  restricted  to 
the  immediate  vicinity  of  the  failed  area  while  tensile 
thermal  stresses  are  distributed  o/er  most  of  the  unfailed 
plate  area. 

Slightly  less  residual  strength  reduction  occurs  if 
the  structurally  failed  area  is  assumed  to  be  immediately 
removed,  forming  a hole,  than  for  the  assumption  that  the 
failed  area  is  retained  .and  continues  to  absorb  and  conduct 
heat.  If  the  structurally  failed  area  is  retained,  subse- 
quent melting  produces  a hole.  Melting  does  not  signifi- 
cantly affect  residual  tensile  strength  except  for  large 
heat  flux  densities  which  produce  melting  at  short  heating 
times.  In  this  case  the  rapid  hole  formation  retards  re- 
sidual strength  reduction  because  most  of  the  subsequent 
incident  flux  passes  through  the  melted  hole. 

The  effect  of  increasing  beam  diameter  is  to  retard 
structural  failure  and  melting,  and  to  increase  residual 
strength  reduction  at  the  longer  heating  times.  Increasing 
beam  power  and.discreasihg  target  plate  thickness  have  the 
’direct  effects  of  promoting  earlier  structural  failure  and 
melting,  and  increasing  residual  strength  reduction. 


The  most  significant  difference  in  the  damage  associat- 


ed with  the  limiting  one-dimensional  heat  flux  cases  is  in 
the  size  of  the  damage  area.  Radial  flux  heating  has  the 
potential  for  producing  failed  areas  larger  than  the  inci- 
dent beam  spot  size  while  axial  flux  damage  is  essentially 
restricted  to  the  area  under  the  beam.  Additional  damage 
due  to  thermal  stresses  outside  any  failed  area  is  produced 
by  radial  flux  in  contrast  to  axial  flux  heating. 

Recommendations 

Based  on  the  results  of  this  study,  the  follov;ing 
recommendations  are  made. 

1.  Additional  NASTRAN  analyses  should  be  conducted  to 
verify  the  one-dimensional  heat  conduction  limits  based  on 
dimensionless  power  per  unit  thickness,  P^g. 

2.  The  apparent  correlation  of  damage  parameters  v;ith 

should  be  more  firmly  established  by  additional  studies  , 

using  the  analytical  tecnniques  and  numerical  radial  code  1 

I 

developc.'d  in  this  study.  i 

i 

3.  The  realism  and  accuracy  of  the  TSTRESS  radial  i 

flux  code  should  be  improved  by  the  following  improvements.  i 

a.  Include  temperature  dependent  material  proper- 

I 

ties.  i 

! 

b.  More  realistic  modeling  of  the  thermal  degrada- 

. tion  of  structural  strength  properties  based  on  laser  i 


heating  exposure  data* 


V 


. an4- convection  i 

d;  XhoXude  Bielt  dynamici  end  interaction  with  the 

* « X 

incident  flux. 

' • * 

4..  Additional  .^terials  should  be  studied. 

5.  iiihear  elastic  fracture  mechanics  should  be  applied 
to  i^roVe  the  realism  of  the  damage  tolerance  analyses. 

6.  The  effect  of  target  plate  size  should  be  analyzed. 

^ 7.  Plasticity  and  target  pre> loading  eiFfects  should  be 

studied..  - - • ' 

' 4...  The  effects  of  assuming^  uncoupled  quasi-static 
heat  conductioh/thermoelasticity  theory  should  be  inves- 
tigated.  ^ 
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numerical  Integration  Algorithm 

Numerical  integration  is  applied  in  subroutines 
STRESSR,  STRONG,  end  RESTRNG  to  c(»npute  thermal  stresses 
and  various  residual  strengths.  In  all  cases,  the  trape* 
zoidal  rule  is  used  (Ref  24:229),  which  gives 

/^ydx  - I 


where 


a,b  « Integration  limits 
• y * f (x) 

h K segment  length  used 

a < X < b 
«■» 

n « (b-a)/h  * number  of 
val  is  divided 


to  divide  the  interval, 
segments  into  %/hich  inter 


In  FORTRAN  code,  a typical  integration  routine  would  be  as 
follows. 


SUM  s 0. 

DO  50  I * J,M 
50  SUM  * SOM  + 2. *1(1) 

SUMM  - SUM  + Y(J-l)  + Y(M+1) 

AB  ^ SUMM  * H/2. 

Numerical  Computation  of  Thermal  Stresses  near  r <■  0. 

For  a solid  plate,  the  thei-mal  stresses  are  given  by 
egns  ($7)  and  (58).  At  the  plate  center  where  r •>  0,  these 
eqns  become  indeterminant.  Hence « separate  algorithms  are 

Us 


necessary  to  compute  thermal  stresses  near  r « 0.  Assuming 
finite  tenqperatures,  an4  applying  I'Hopital's  rule  to  eqns 
(57)  and  (58)  gives  the  following. 

* 006  (0)  * oE  l.r/**Trdr-i^T(0)j  (110) 

This  equation  is  numerically  integrated  in  subroutine 
STRESSR  to  obtain  stresses  at  r « 0. 

It  was  also  found  that  a special  algorithm  was  required 
for  computing  stresses  at  the  first  radial  increment  away 
from  r » 0 because  eqns  (57)  and  (58)  are  still  erroneous 
for  small  values  of  r.  The  following  texm  in  these  equa- 
tions becomes  erroneously  large  for  very  small  values  of  r. 

r 

i.,  / Trdr 
r*  0 

To  evalupce  this  integral  at  the  first  radial  increment, 
repl-^  -e  the  upper  limit  with  r « AR,  where  AR  = radial  ele- 
Kk'.at  width.  Assuming  a linear  temperature  distribution  be- 
tween r « 0 and  r » AR, 

T(r)  - T(0)  [l-Pr]  (111) 


where 

T(0)  temperature  at  r ■>  0 

P «•  constant  (slope  of  temperature  distribution) 


116 


Define 


.AR  AR 

ARR  - / rTdr  « T(0)  / r(l-Pr)dr 
0 0 


(112) 


The  value  of  P is  foun^  as  follows. 


T(l)  « T(0)  [l-P(AR)] 


(113) 


where 


T(l)  * temperature  at  r « AR 


P » 


SI  TOTJ 


(114) 


o 


Substituting  into  eqn  (112)  and  integrating  gives 


ARR  » (AR)Mt(0)/6+T(1)/3J 


(115) 


This  algorithm  is  used  in  computing  thermal  stresses  at  the 
outer  radius  of  the  first  radial  element  for  a solid  disk. 


O 


Radius  of  Melted  Hole  (RMELT)^ 

> 

An  element  is  considered  to  be  melted  (transformed 
fr^  the  solid  to  the  liquid  ^ase)  when  the  temperature  of 
. that  element  reaches  the  upper  melting  temperature  (TUMELT) 

^ which  is  638*C  for  2024-T3  aluminum.  The  phase  change  is 
accounted  for  as.  described  in  Section  VI.  The  assumption 
is  iftade  that  the  melt  is  removed  instantaneously  upon  reach- 

u,  ■ ' 


ing  the  upper  melting  temperature.  Since  the  axially  sym- 
metric beam  is  considered  to  be  located  over  the  center  of 
a flat  plate,  melting  produces  a circular  hole  in  the  plate 
which  progressively  grows  radially  with  continued  applica- 
tion of  heat  flux.  The  hole  growth  ceases  when  the  hole  be- 
comes so  large  that  any  remaining  irradiated  plate  area  ab- 
sorbs insufficient  heat  flux  to  produce  further  melting. 

The  radius  of  the  melted  hole  is  given  by 

RMELT  I*DELR  (116) 


where 


I a element  number  of  outermost  element  which  has 
been  heated  to  upper  melting  temperature  (638**C) 

DELR  a element  radial  width 

i 

Radius  of  Failed  Area  due  to  Thermal  Degradation  (RFD) 

This  quantity  is  computed  in  Subroutine  RADPAIL  (lines 
18  to  21).  Jt ‘is  postulated  that  a failed  area,  centered 
at  the  coincident  beam  and  plate  cehters,  is  generated  or 
expanded  whenever  any  of  the  three  structural  strength  pro- 
perties is  reduced  to  zero,  as  defined  by  Figures  2 through 
* ' - ■* 

4.  As  noted  in  Section  VI,  all  three  strength  properties 
of  an  element  are  talcen  to  be  zero  if  the  temperature  of 
that  element  is  greeted  than  370 *C  (for  2024-T3  aluminum). 
Hende,  the  radius  of  a failed  area  due  to  thermal  degrada- 
tion  is  considered  to  be  dependent  on  temperature  only  and 
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ia  developed  in  a manner  analagous  to  that  of  a hole  due  to 
melting  (RMELT) . That  is,  the  area  radius  given  by 


RPD  » I*DELR  (117) 


inhere 

X * element  number  of  outermost  element  which  has 
been  heated  to  zero  structural  strength  temp- 
erature (370*C) 

Radius  of  Failed  Area  due  to  Thermal  gtresses  (RFS) 

This  quantity  is  computed  in  Subroutine  RADFAlli  (lines 
22  to  34)  and.  is  based  on  the  postulate  that  a beam  axis 
^centered  failed  area  is  generated  in  tlxe  heated  plate  if 
^ j any  of  the  three  thermal  strels  components  exceeds  its  re- 

spective room  temperature  strength  criterion.  This  para-* 

meter  is  used  as  a measure  of  the  failed  area  that  would  be 

* 

produced  by  thermal  stresses  independent  ot  the  action  of 

■»  a 

structural  strength  thermal  de^jradat^on.  The  thermal  stres- 
ses of  each  element  are  compared  to  room  temperature  ulti- 
mate  tensile  strength*,  compressive  yield  strength,  and  ulti- 
mate shear  strength.  If  any  of  ttwise  criteria  are  exceeded 
at  any  element,  the  plate  material  is  considered  to  have 

. t 

zero  strength  at  the  element.  A failed  area  radius  Is  then 
defined  by  the  outermost  radial  element  which  has  reached  a 
failure  strength  level*  ^ ^ 

O • 
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Radius  of  Failed  Area  due  to  Combined  Action  of  Thermal 
Stresses  and  Degradation  (rfds) 

This  quantity  is  conqputed  in  Subroutine  rAdpail  (lines 

35  to  70)  and  is  based  on  the  postulate  that  a radial  ele> 

sent  is  reduced  to  2e.r0  strength  if  any  of  the  three  stress 

conponents  exceeds  its  respective  elevated  temperature 

structural  strength  as  defined  in  Section  VI.  The  radius 

of  a failed  area  due  to  the  combined  action  of  thermal 

stresses  and  thermal  degradation  of  structural  strength  is 

determined  by  the  outer  radius  of  the  outermost  element 

which- has  one  or  more  stress  components  greated  than  its 

failure  strength. 

Actual  Radius  of  Failed  Area  (RFAIL) 

Two  options'  cohcerniiig  plate  failure  zones  are  provided 
in*^Program  TSTRESS,  Although  structurally  failed  areas  are 
defined  by  RPSr  RFD,  and' RFDS,  it  Us  cohceivable  that  these 
failed  circular  areas  might  either  remain  in  place  in  the 
plate  or  be  immediately  removed  by  some  source  such  as  aero- 
dynamic pressure.  In  either  case  it  is  postulated  that  any 
melted  zone  (RMELT)’  is  immediately  removed.  Hence,  one 

option  considered  (defined  by  IR  * l,'^  line  108  of  TSTRESS) 

► ’ 

is  that  of  allowing  a hole  to  be  formed  by  either  of  the 
mechanisms 'defined  jby  RNSLT  and  RFDS.  That  is,  th%  actual 
h61e  radius  (RFAIL)  would  be  the  larger  of  the  above  radii, 

ln„  which  ease  immediate  removal  of  this  zone  is  assumed  re- 

^ ' *"  * 

gardless  of  f)ie  failure  ipechanismi.  The  second  option  (de- 
fined  by  IR  » 2)  removes  the  failed  area  only  if  the  failure 


nechanisn  is  Belting.  I^t  is  the  actual  hole  radius  is  al> 
ways  equal  to  RMELT,  in  which  case  laser  beeun  heat  flux  con* 
tinues  to  be  absorbed  by  the  plate  area  within  a structural- 
ly failed  circle  and  outside  the  melted  hole.  As  shown  in 
the  sample  output  of  TSTRESS,  Appendix  A,  for  either  of  the 
2dt>ove  options,  all  of  the  «ux>ve  failed  r^dii  are  listed. 
Following  the  failed  radii  listing  are  listings  of  the 
stress  ccxnponents  which  caused  failure  due  to  thermal  stres- 
ses alone,  and  due  to  combined  thennal  stresses/degradation. 

Failed  Hole  due  to  Stress  Concentration  (RFSC) 

Following  the  above  listings  in  the  sample  output,  Ap- 
pendix  A,  is  a table  of  failed  hole  radius  due  to  stress 
concentration  (RFSC)  as  a function  of  the  remotely  applied 
tensile  stress  (SIGAPP)  in  the  plane  of  the  plate.  The 
initial  hole  radius  is  taken  as  that  due  to  structural  fail- 
ure, It  is  assumed  that  failure  due  to  melting  is  associat- 
ed with  high  temperatures  at  the  edge  of  the  hole  which  tend 
to  soften  the  material  to  the  extant  that  stress  concentra- 
tions are  negligible.  For  a plate  in  uniaxial  tension,  the 
r^tio  of  local  tersile  stress  to  the  remotely  applied  ten- 
sile stress  is  (Ref  21:399) 

* 1 + in  + i.  n*  ai8) 

C • t t 

•• 

where  ’ s,  • , 
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n * a*A*^. 


a B hole  radius  ^ RPSC 

j * ronotely  applied  tensile  stress  » SIGAPP 

00  « local  tensile  stress  9r 

If  0*0  exceeds  roon  temperature  ultimate  tensile  strength, 
the  failed  radius  (RPSC)  is  computed.  Computations  are 
made  foi  several  values  of  remotely  applied  stress  (SIGAPP) . 
These  calculations  are  made  in  Subroutine  RADSC  amd  in 
TSTRESS  at  lines  238  to  240. 

Residual  Tensile  Strength  Due  to  Thermal  Degradation  (RESDEG) 
The  uniformly  distributed,  remotely  applied  tensile 
force  required  to  fail  a solid  plate  in  tension,  consider" 
Ing  only  the  effects  of  thermal  degradation  of  ultimate 
^ tensile  strength  is  given  by  (Ref  31) 

^^0(t)'l-6dx  ~ ^ {OQ+M{T(x,t)-ToI  }dzdx  (119) 

where 

o(t)  « STRNGTH  «=  RESDEG  = residual  tensile  strength 
due  to  thermal  degradation. 

Oq  **  ■»  original  unheatev!  tensile  strength 

M » 'Slope  of  thermal  degradation  plot  of  ulti-* 
n^if  tensile  strength  vs.  temperature 

T(x,t)'  » Ibcal  temperature  of  pic^te  at  dxstance  x from 
plate  center 

T0  w temperature  corresponding  to 
2b  » • plate  width 

0 ^ 1 * plate  thic)^ess  ' 
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Perfonaing  the  indicated  integrations  on  the  left  side  and 
integrating  with  respect  to  z on  the  right  side  gives 


1 

o(t)  * ~ / {Oo+M[T(x,t)-Tj}  dx 

Due  to  axial  synanetry,  this  is  equivalent  to 

, b 

o(t’  “ g f {aQ+MCT(x,t)-Tp)}  dx 


(120) 


(121) 


As  discussed  in  Section  VI,  the  thermal  degradation  curve 
was  approxiisated  by  two  straight  lines.  Eqn  (121)  then  be- 
comes 

; 

b* 

o(t)  * I I ^ {0Q2+M2l.T(x,t)-T<52J>  dx 

b ; , \ 

+/  (T(x,t)-ToiJ)  dx>  (122) 

b*  ' 

where  b’*’  = value  of  x where  T(b*,t)  » Tqz 

Toy  « temperature  at  inter  .-section  of  the  two  straight 
' line  approximatxone 

My  ■ slope' of  degradation  cufve  for  370®C  > T(x,t)  > 

7*02 

« 

Ml  * slope  of  dsgradiition  curve  for  < T(x,t)  < 

' , Toy 

Toi  - To 

®oi  " ^6  “ ^to 

002  •*  ultimate  tensile  strength  at  To2, 

Allowance  a failed  area  ar  cae  pl^te  center  is  accounted 


for  by  beginning  the  integration  at  the  edge  of  this  area. 


o(t) 


1 

b 


b* 

./  IOo2+MiT{x,c>~MjTo2ldx 


b 

/ IOoi+MiT(x,t)-MiTo,) 

« A a 


idiere  bf  « BF  *■  raoius  of  failed  area 
The  following  approximations  also  are  made: 


(123) 


Oq,  at  0.75  Oqj 
To»  » 0*C 


This  results  in  the  tom  of  the  algorithm  used  in  subroutine 
STRONG. 


o(t)  » i|(b-0.25fc*-0.75bf)Oj„ 
b ' 


b* 

-Ma (Tq2 (b*-bf)-  / T(x#t)dx) 
bf 

b 

+Mi  / T(x>t)dx 

b* 


w^re’  Mj;  » -200  for  2024-T3  aluminum 

Mj  -'-55 

<^0i  - Ho  - .«t820  n/cm* 

T^jj  « 20O*C 


(124) 


the  proportion  of  the  original  urihea ted  ultimate  tensile 


strength  is  then  computed  as 


POBG  « P • STRNGTH/Ojjj  (125) 

Residual  Tensile  Strength  Due  to  Melting  (RESMELT) 

This  quantity  is  a direct  function  of  the  reduction  in 
plate  cross-sectional  area  due  to  melting  of  a hole,  and  is 
computed  in  the  main  program  (lines  245-246)  as  follows 

PMELT  » (1-RMELT/B)  (126) 

RESMELT  « PM£LT*ULTENS  (127) 

where 

PMELT  « proportion  of  original  unheated  ultimate 
tensile  strength  remaining 

RMELT  “ radius  of  melted  hole 

B plate  radius 

ULTENS  » = original  unheated  ultimate  tensile 

strength 

RESKELT  *■  residual  tensile  strength  due  to  effect 
of  melted  hole  only 

Residual  Strength  Due  to  Thermal  Stresses  (RESTRSS) 

The  tensile  load  carrying  capacity  of  the  rectangular 
platp  at  any  time  t,  considering  only  the  effects  of  thermal 
Stresses  and  a failed  area  at  the  plate  center,  is 


P » P^u  - ®a  t2Jl(b-b«))  - Ptu 


h. 

b 


(128) 


where 


* rooM  ten^rature  ultimate  tensile  load  strength 

■ thermal  tensile  stress  averaged  over  area  o£ 
applied  load 

t *B  plate  thickness 

b • plate  radius  or  width/2 

* radius  of  failed  area 

bs  « radius  at  which  tensile  stress  (09)  becomes 
greater  than  zero 

b«  * the  greater  of  b«  or  bt 

Dividing  by  the  original  plate  cross-sectional  area,  2bt, 
gives  the  residual  tensile  strength  due  to  thermal  stresses 
only. 


(129) 


where 


^to  “ ULTENS  « room  temperature  ultimate  tensile 
strength 

0 **  RES7RSS  * residual  tensile  strength  due  to 
thermal  stresses  only 


Dividing  by*  gives  the  proportion  of  the  original  tensile 
strength  remaining. 

PSTRSS  « RESTRSS/OLTENS 


The  average  tensile  stress  is  computed^ by  numerical  inte- 
gration of 


o 


(130) 


1 l> 

* AVGfTEMS  • O0df 


Residual  Tensile  Strength  Due  to  the  Combined  Actions  of 
tieltinq»  Thermal  Stresses » and  Thermal  Degradation  (RESTRS2) 

* 

This  algorithm  is  similar  to  that  for  thermal  stresses 
except  that  the  elevated  temperature  ultimata  tensile 
strength,  averaged  over  the  net  section  area,  replaces  the 
room  temperature  ultimate  strength. 

o - ?tu  - 0,  (131) 


tdiere 


0 » RSSTRS2  « residual  tensile  strength  due  to 

melting,  therisal  stresses  and  thermal  degrada- 
tion 

Ftu  “ ULTEMP  » elevated  temperature  ultimate  tensile 
strength  averaged  over  net  section  area 

hr  e « B7  > radius  of  failed  area 


ULTEMP  is  computed  by  numerical  integration  of  the  thermal 
degradation  equations  from  Section  VX  for  ultimate  tensile 
stxviihgth. 


ULTEMP  - 


b 

^^Ftudr 


(132) 


i 

where 

■ < ' 

» • 

<B  ULT(I)  w ^lev^ed  temperature  ultimate 
tensile’  strength  . 

, '>■ 

The  proportion  of  room  tes^erature  ultimate  strength  is 
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PSTRS2  • 8ESTRS2/DLTENS  (133) 

o 

Hesiduyl  Tensile  Strength  Due  to  Failed  Area  (RESPAIL) 

The  above  neasures  of  residtaal  strengths  contain  the 
effects  of  a reduction  in  net  section  load-carrying  area 
due  to  melting  or  structural  failure.  The  effect  of  the 
structural)/  failed  area  on  residual  tensile  strength  is 
coB^uted  ija  the  main  program  (lines  249-250)  as  follows. 

PO  « (1-RPAIL/R)  (134) 

RESPAIL  » IO*ULTENS  (135) 

where 

(..) 

PO  « proportion  of  room  temperature  ultimate  tensile 
strength  remaining 


O 
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MASTRAN  Validation  of  Prograitt  TSTRESS 

To  validate  the  one-dinensional  radial  heat  conduction 
program,  teiq>erature  and  thermal  stress  distributions  were 
ccmapared  with  those  frcsa  a two~dimensiOnal  finite  element 
analysis  using  HASTRAN  (Ref  2) . The  MASTRAN  finite  element 
model  consisted  of  a segment  of  a solid  cylinder  as  shown 
in  Fig.  28.  The  wedge  shaped  segment  was  divided  into  15. 
radial  elements  and  5 depth  elements.  The  radial  width  of 
the  elements  varied  from  a fine  mesh  under  the  laser  beam 

I 

to  a relatively  broad  mesh  at  the  outer  radii.  The  element 
thicknesses  were  constant,  each  being  1/5  of  the  plate 
thickness. 

The  fledge  segment  is  sufficient  to  model  the  problem 
of  interest  because  of  the  axial  symmetry  of  the  heat  flux 
input  which  implies  that  the  problem  is  independent  of  the 
radial  coordinate.  The  wedge  angle  was  chosen  to  he  15 
degrees  to  provide  reasoneble  aspect  ratios  (width/length) 
for  the  finite  elements.  The  radius  of  the  wedge  was  taken 
as  15  cm  to  be  consistent  with  the  one-dimensional  numeri- 
cal radial  flux  model.  In  general,  this  study  was  restrict- 
ed to  large  plate  diameter/beam  diameter  ratios  to  minimize 
the  effect  of  this  parameter  on  the  results.  Also,  as  in 
the  numerical  radial  model,  no  heat  losses  were  considered 
in  the  MASTRAN  analysis. 

The  heat  flux  input  to  the  surface  elements  under  the 
beam  radius  was  approximated  by  using  the  flux  density  value 
at  the  radial  midpoint  of  each  surface  element.  For  a 


Pig.  28  Two-Dimensional  Finite  Element  Model 


Gaussian  beeun,  the  heat  flux  density  distribution  is  given 
by  eqn  (9)  as 

-2rVa* 

I(r)  « Ipae 


where 


a « 2o  *>  bean  radius 

This  is  the  sane  approximation  used  for  the  heat  flux  input 
to  the  one->dimensional  model.  However,  the  one"dimensional 
nodel  has  constant  width  (radial)  elements  sized  to  achieve 
a cmi^romise  between  accuracy  and  computational  efficiency. 
This  was  initially  accomplished  by  analyzing  the  effect  of 
reducing  element  width  on  computer  time,  and  on  temperature 
and  thermal  stress  values.  An  element  width  of  0.3  cm  was 
selected  as  a reasonable  compromise  based  on  the  judgment 
that  larger  widths  gave  significantly  different  temperatures 
and  thermal  stresses  while  smaller  widths  produced  very 
small  differences  in  temperatures  and  thermal  stresses  but 
rapidly  increased  computational  cost. 

Figures  29  and  30  show  cos^arisons  of  thermal  stresses 
between  the  one-dimensional  radial  model  and  the  NASTRAN 
two-dimensional  finite  element  results.  Such  comparisons 
were  made  for  four  different  input  conditions*  involving  two 
different  heat  flux  densities  and  two  plate  thicknesses. 

Only  one  beam  diameter  was  considered,  3 cm..  The  results 
presented  in  Figures  29  and  30  are  typical.  The  NASTRAN 


RADIAL  STRESS  (W/CM*)xlO 


stresses  are  at  the  heated » or  fronts  surface.  For  the  con-* 
ditlons  listed  in  the  figures#  the  heat  conduction  is  pri- 
Mrily  radial  as  determined  by  comparing  front  and  back  sur- 
face temperatures  and  stresses  as  partially  listed  in  Table 
I.  The  relatively  small  front-to-back  differences  in  the 
temperatures  indicate  that  the  axial  heat  flux  vector  com- 
ponent is  small  since  the  heat  conduction  is  directly  pro- 
portional to  the  spatial  thermal  gradient  (reference  eqn  (21)). 

Two  conclusions  drawn  from  Fig.  29  are  that  the  one- 
dimensional radial  flux  model  does  not  give  zero  slopes  in 
the  stress  distributions  at  r » 0#  and  that  the  radial  flux 
model  stresses  appear  to  increase  (become  more  negative) 
relative  to  the  NASTRAN  stresses  with  increasing  time. 

These  trends  also  are  exhibited  in  Pig.  30.  Zero  slope  at 
r * 0 is  required  by  the  axial  symmetry  of  the  problem. 

In  order  to  evaluate  the  significance  of  these  trends# 
the  temperature  distributions  were  compared  in  greater  de- 
tail as  shown  in  Fig  31.  Solid  curves  were  fitted  through 
the  NASTRAN  data#  and  the  one-dxmensional  radial  flux  model 
data  for  the  same  conditions  arv^  represented  by  the  triangle 
symbol.  The  same  trends  exist  in  the  temperature  compari- 
sons as  noted  for  the  previous  stress  comparison?. 

i The  expanded  radius  scale  allows  all  of  the  NASTRAN 
grid  point  values  to  be  plotted.  In  Figures  29  and  30 
several  of  the  grid  point  values  near  r «•  0 were  not  plotted 
in  order  to  clarify  the  plot.  Since  the  temperature  compar- 
isons showed  similar  results  to  the  stress  comparisons#  the 
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Table  I • NASTRAN  Fjfonrt  and  Back  Surraca 

Temperatures  and  Stresses  at  Selected 
Radii  (for  conditions  listed  in  Pig. 


possibility  that  the  differences  noted  are  due  to  the  dif- 
ferences in  element  size  was  investigated.  As  shown  by  the 
data  in  Pig.  31,  the  resolution  of  the  NASTRAN  finite  ele- 
ment model  was  much  finer  near  the  beam  axis  (r=0)  than  that 
of  the  constant  width  elei..ent  raoial  flux  model.  For  the 
data  presented  thus  far,  the  ccxnputational  time  increment 
was  0.05  sec.  in  both  clie  NASTRAN  analysis  and  the  one- 
dimensional radial  flux  results.  In  order  to  reduce  the 
element  size  in  the  one-dimensional  radial  model  to  that  of 
the  NASTRAN  model  under  the  beam,  the  heat  conduction  stabil- 
ity criterion  of  eqn  (28)  requires  a consequent  reduction  in 
computational  time  increment.  These  combined. reductions  re- 
sult in  a very  large  increase  in  computer  time  cost  (on  the 
order  of  a factor  of  20) . 

The  element  size  for  the  one-dimensional  radial  model 
was  reduced  to  0.05  cm,  requiring  that  the  time  increment  be 
reduced  to  0.00^  sec,  and  the  results  are  represented  in 
Fig.  31  by  the  square  symbols.  Comparison  with  the  previous 
radial  model  data  in  the  seime  figure  shows  that  the  smaller 
element  size  does  develop  a flatter  slope  near  r » 0.  Hence, 
it  is  concluded  that  the  apparent  lac)c  of  zero  slopes  in  the 
one-dimensional  radial  modal  data  is  due  to  the  relatively 
coarse  spatial  resolution  used  to  achieve  a compromise  with 
computer  run  time. 

The  trend  for  the  one-dimensional  radial  model  value.3 
to  increase  with  time  relative  to  the  NASTRAN  values  is 
still  evident.  The  one-dimensional  radial  model  with  in- 
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creased  resolution  is  seen  to  give  slightly  reduced  values 
ar  r B 0 while  those  at  the  larger  radii  are  increased  es> 
specially  for  the  longer  heating  time.  This  result  is  be- 
lieved to  be  least  partially  attributable  to  plate  dia- 
meter effect  . he  one-dimensional  radial  model  code  was 
developed  with  a one-dimensional  array  storage  limit  of  100 
spaces  which  was  judged  to  be  a reasonable  upper  limit  on 
the  number  of  radial  elements  to  be  applied.  For  the  anal- 
ysis plotted  in  Fig.  31  with  element  widths  of  0.05  cm,  this 
limited  the  plate  radius  to  5 cm.  Recalling  that  the 
NASTRAN  data  are  for  a plate  radius  of  15  cm,  it  is  reasoned 
that  the  effect  of. the  insulated  (no  heat  losses)  plate  edge 
is  being  felt  and  causing  temperatures  to  increase  at  the 
larger  radii.  If  this  is  the  case  it  would  be  expected  that 
this  effect  would  cause  the  outer  radii  temperatures  rela- 
tive to  those  from  KASTriA.N  to  increase  with  time.  It  would 
be  expected  that  this  effect  would  propagate  to  smaller 
radii  with  incireasing  time. 

The  data  in  Fig.  31  are  consistent  with  the  expected 
effects  of  reducing  the  plate/beam  dieuneter  ratio,  although 
it  is  not  certain  whether  this  effect  is  significant  enough 
to  be  causing  the  one-dimensional  radial  5 second  heating 
time  values  to  exceed  the  NASTR/iN  values  near  r * 0. 

Fig.  32  presents  additional  comparisons  bf  temperature 
distributions  for  a plate  thickness  of  0.3  cm.  In  this  case 
the  heat  flux  vector  contains  a slightly  larger  axial  compo- 
nent than  for  the  0.12  cm  thick  plate  at  corresponding  heat- 
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ing  times.  This  is  indicated  by  the  comparison  of  the  front 
and  back  surface  teng>erature8  from  the  NASTRAN  analysis  as 
shoi«n  in  Fig.  32. 

Comparison  of  the  StASTRAN  front  surface  and  radial  flux 
temperatures  shows  the  same  trends  as  noted  in  Fig.  31  for 
the  thinner  plate.  The  effect  of  reducing  the  radial  model 
element  width  from  0.3  cm  to  0.05  cm  also  was  evaluated  for 
the  thicker  plate  as  shown  in  Fig.  33.  Again  the  results 
are  similar  to  those  of  Fig.  31  and  no  additional  insight 

is  provided  on  the  possible  cause  of  the  increase  with  time 

« 

of  the  radial  model  ten^eraturcs  near  r « 0 relative  to 
those  from  the  NASTRAN  analysis. 

Based  on  the  preceding  comparisons,  it  is  concluded 
that  the  one*’dimensional  radial  flux  munerical  model  pro* 
duces  temperature  distributions  which  are  generally  within 
10  percent  oi  those  from  the  two-dimensional  NASTRAN  analy- 
sis. Peak  compressive  stresses  also  are  in  close  agreement 
between  the  two  analyses.  Tensile  stress  distributions  are 
in  reasonable  agreement  although  conclusions  regarding  peak 
tensile  stresses  are  not  definite  due  to  inadequate  resolu- 
tion in  the  NASTRAN  finite  element  model  at  the  larger  radii 
(see  Fig.  30) . Improved  evaluation  of  the  radial  flux  model 
could  be  obtained  by  modifying  the  radial  flux  code  to  in- 
crease the  maximum  number  of  elements  from  100  to  300  which 
would  eliminate  the  plate  size  effect  introduced  by  decreas- 
ing the  element  width  to  0.05  cm.  Also  the  NASTRAN  finite 
element  model  (Fig.  28)  should  have  additional  radial  ele- 
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33  Temperature  Distribution  Comparisoni  Effect  of  1-D 
Radial  Heat  Flux  Model  Element  Size 
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■ents  added  to  ioqprove  the  resolution  in  the  areas  where 
peak  tensile  stresses  occur;  and  NASTRAN  data  should  be 
geneiated  for  longer  heating  tians  than  5 seconds  to  better 
evaluate  the  possible  trend  for  the  radial  model  tempera- 
tures to  increase  with  time  relative  to  the  NASTRAN  tempera 
tures. 
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